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Abstract 

A short review of all available results (perturbative, nonperturbative, and 
exact) on (i-dimensional spin models is presented in order to introduce the dis- 
cussion of their 1/N expansion at d = 2, where the models are asymptotically 
free. 

A general two-dimensional spin model with U(iV) invariance, interpolating 
between CP^ -1 and 0(2N) models, is studied in detail in order to illustrate 
both the general features of the 1/N expansion on the lattice and the specific 
techniques devised to extract scaling (field-theoretical) behavior. 

The continuum version of the model is carefully analyzed deriving quan- 
titative 0(1/N) physical predictions in order to establish a benchmark for 
lattice computations. 

The 1/N expansion on the lattice, including second-nearest-neighbor in- 
teractions, is set up by constructing explicitly effective propagators and ver- 
tices, and exhibiting a number of exact results and integral representations 
that allow a substantial reduction of the numerical effort. The technique of 
asymptotic expansion of the lattice propagators, basic to the derivation of 
analytical results in the scaling domain, is presented in full detail and applied 
to the model. Physical quantities, like the free energy and different definitions 
of correlation length, are evaluated. 

The lattice renormalization-group trajectories are identified and univer- 
sality among different lattice (and continuum) schemes in the scaling region is 
explicitly proven. As a byproduct, representations of the 0(1/N) contribution 
to the yl-parameter ratios and to the lattice (3 functions are obtained. 

A review of other developments based on the lattice 1 /N expansion (finite 
size scaling, higher orders, fermionic models) is presented. 



PACS numbers: 11.15 Ha, 11.15 Pg, 75.10 Hk. 

To appear in "Rivista del Nuovo Cimento' 
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... this (is) the statement of the research, lest human deeds should in 
time fade, and great and wonderful works... become unknown... 



[Herodotus, Histories, 1, 1] 



A better qualitative and quantitative understanding of quantum field theories requires 
an improvement of the analytical and numerical methods of approximation. Lattice field 
theories are a natural ground of application of large-scale numerical techniques. More effi- 
cient algorithms and more powerful computing machines lead to an ever increasing amount 
of numerical results. These results are however affected by two major limitations: 

1. the quality of the information grows as the logarithm (or, at best, as a small power) 
of the numerical effort; 

2. the lack of control on the systematic errors possibly induced by some of the numerical 
techniques is often tantalizing. 

Both limitations are intrinsic to the field-theoretical nature of the problem addressed: off- 
critical systems, and even some bulk properties of critical systems, can usually be studied 
with great numerical precision as well as with sensible analytical methods. Systems at 
criticality, and the extraction of their scaling properties, however constitute a much more 
difficult challenge. 

In view of the abovementioned limitations, the parallel development of more powerful 
analytical techniques, even with a limited domain of applicability, is certainly very welcome. 
On one side it allows the comparison between numerical an analytical results in a controlled 
environment, and therefore it leads to a check of applicability for those techniques and 
methods whose reliability cannot be taken for granted a priori. On the other side it may 
improve our understanding of those properties that cannot be directly tested with present- 
day numerical methods and may strengthen some of the theoretical hypotheses that must 
unavoidably be used in the field-theoretical interpretation of numerical data. One can easily 
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understand the need for such theoretical pieces of evidence when one realizes that, notwith- 
standing all perturbative results and the substantial agreement existing in the theoretical 
physics community, there is at present no independent nonperturbative proof of existence of 
any asymptotically-free quantum field theory, whose relevance in the construction of models 
of the fundamental interactions cannot be overstressed. 

In this perspective, we think that the 1/N expansion |3J (expansion in the number 
of field components) may be a rather rewarding instrument of analysis not only in the 
context of continuum field theories, where it has long been known as a major source of 
nonperturbative information, but also in the case of lattice field theories, where to the best 
of our knowledge the 1/N expansion is the only approach leading to some theoretical evidence 
for the existence of a continuum limit and of a nonvanishing scaling region, where the field- 
theoretical properties of the models can in principle be explicitly tested with predictable 
precision. 

The conceptual foundations of the power of the 1/N expansion are essentially the fol- 
lowing: 

1. iV is an intrinsically adimensional parameter, representing a dependence whose origin 
is basically group-theoretical, and leading to well-defined field representation for all 
integer values, hence it is not subject to any kind of renormalization; 

2. iV does not depend on any physical scale of the theory, therefore we may expect it to 
play no role in the parametrization of criticality. As a consequence there is no physical 
reason not to expect reasonable convergence properties from an expansion in 1/N, at 
least in well-defined regions of the other physical parameters. 

Evidence for a finite radius of convergence of the 1/N expansion has been produced 
in a number of instances, notably in the proposed exact S'-matrices for a number of two- 
dimensional bosonic and fermionic models. More generally, the large order behavior of the 
coefficients of the 1/N expansion can be studied by applying inverse-scattering techniques 
to the problem of finding instanton solutions of the effective actions In the case of 

0(A r )-symmetric (^f 2 ^ theories in less than four dimensions, the 1/N perturbation series 
can be resummed by a Borel transformation, and in the two-dimensional nonlinear a model 
one is led to conjecture the convergence of the series also for the Green's functions. 

Till now the major domain of application of the 1/N expansion has been in the evaluation 
of critical exponents for many different classes of models in the range of dimensions com- 
prised between 2 and 4 (lower and upper critical dimensions). At the critical dimension the 
critical exponents are trivial, but the logarithmic deviations from scaling and the dynamical 
mass generation lead to the rich phenomenology characteristic of asymptotic freedom. This 
phenomenology can be studied, on the lattice as well as in the continuum, by applying the 
1/N expansion in conjunction with proper modifications of the methods usually adopted 
in standard renormalizable quantum field theories, notably the techniques of regularization 
and renormalization of the physical parameters. Due to the nonrenormalized character of N, 
the 1/N expansion leads to results whose renormalization-group invariance properties are 
much more transparent than those of standard perturbation theory. However, as far as we 
could check, the 1/N expansion commutes with perturbation theory, and therefore it pro- 
vides a direct reinterpretation and an unambiguous resummation of the perturbative results. 
Such phenomena as the effects of a change of regularization scheme, the role of dimensional 
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transmutation in the parametrization of renormalization-group invariance, the relationship 
between dynamical mass generation and Borel ambiguity in resummations, the interplay be- 
tween dimensional regularization, minimal subtraction, and £-expansion, the mixing effects 
and the subtraction of perturbative tails in the evaluation of quantum expectation values 
of composite operators, all find a specific and transparent illustration in the context of the 
1/N expansion already at the lowest nontrivial order of computations |5|||. 

The main purpose of the present review is to describe the results that may be obtained by 
applying the 1/N expansion to the lattice versions of two-dimensional spin models. Therefore 
we shall only briefly sketch the main results presented in the (wide) literature on all the 
abovementioned topics, focusing only on those continuum results that are essential in order 
to introduce their lattice counterparts. 

The subject of lattice 1/N expansion has not till now received a systematic treatment: 
as a consequence many sections of the present paper (especially Sect. § and Sects. |7Hl4[) are 
essentially original work by the authors. 

In order to present a few peculiar techniques and results of the application of the 1/N 
expansion to renormalizable lattice field theories, we found it proper to focus on a specific 
class of models. These models should be simple enough to make all calculations as short 
and understandable as possible, as well as to make accurate numerical simulations feasible 
now or in the near future. However it was necessary for the completeness of the presentation 
to deal with a sufficiently rich phenomenology, notably nontrivial mass spectra, gauge and 
topological properties, besides the obvious request of perturbative asymptotic freedom. After 
much thought, we decided to study a two-parameter model of two-dimensional spin fields 
with U(iV) global symmetry [0]. This model interpolates between the standard 0(2N) 
|8HT0| and CP^ -1 models |11,T^|, and for special values of the parameters it represents the 



gauge-fixed bosonized version of the minimal coupling of massless fermions to CP fields 

HI- 

The present paper is organized as follows: 

In Sect. |2] we briefly discuss the general results that have been obtained in the study of 
(i-dimensional spin models, with special emphasis on the topics related to perturbative and 
nonperturbative renormalizability. 

In Sect. |H we introduce a class of 1/iV-expandable two-dimensional spin models, and 
discuss the qualitative picture of their properties that one may draw from a large- N analysis. 

In Sect. f| we review a number of exact results, especially factorized S'-matrices, that 
apply to the models under consideration for peculiar values of the parameters. 

In Sect. [5] we discuss the 1/N expansion in the continuum version of the models, introduc- 
ing our regularization and renormalization procedure, defining observables, and extracting 
some quantitative 0(1/ N) physical predictions. 

Sect. H is dedicated to the presentation of a number of alternative lattice formulations 
and motivates our choice of a lattice action, which will turn out to depend explicitly on an 
Abelian vector field and on an extra parameter eventually allowing for a Symanzik tree-level 
improvement of the action. 

In Sect. |7] the basic ingredients of the 1/N expansion on the lattice, effective action, 
propagators, and vertices, are introduced. 

Sect. |8] is devoted to a specific technical problem, the search for integral representations 
of the effective lattice propagators that, in the case of nearest-neighbor interactions, allows 
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a substantial simplification of the numerical tasks in the evaluation of the effective Feynman 
diagrams. 

In Sect. |9] we introduce the basic ingredient of all lattice computations in the scaling (field- 
theoretical) regime, the asymptotic expansion of the lattice propagators for small values of 
the (dynamically generated) mass gap, i.e. for values of the correlation length much bigger 
than the lattice spacing (we are always working in the infinite- volume limit). 

In Sect. [It] special techniques for the asymptotic expansion in the case of nearest-neighbor 
interactions discussed in Sect. || are presented. 

Sect. |TT]is dedicated to applying the abovementioned results to the actual 0(1 /N) eval- 
uation of physical quantities in the scaling region of our class of models. We discuss the 
possible definitions of correlation length, show the universality of the lattice results, and 
give a full analysis of the simplest correlation function, the two-point correlator of the fun- 
damental fields, including the evaluation of the lattice wavefunction renormalization factor. 

In Sect. [T2] the abovementioned lattice contributions to physical quantities are explicitly 
evaluated in the case of nearest-neighbor interactions. 

Sect. [L3| is devoted to the issue of topological operators on the lattice; different definitions 
are analyzed and compared in the context of the 1/N expansion of CP*' 1 models. 

In Sect. ^ we rephrase our results in the language of standard perturbation theory and 
perturbative renormalization group. We discuss the evaluation of the ratio of A parameters 
in the context of the 1/N expansion and extract an explicit representation of the 0(1/ N) 
contributions to the lattice renormalization-group (3 function, clarifying some subtleties con- 
cerning the noncommutativity of some limits at the border of the space of parameters, which 
however does not affect the physical predictivity of the model. 

In Sect. 15 we review some results concerning the possibility of performing a finite-size- 
scaling analysis of spin models by the help of 1/iV-expansion techniques. 

In Sect. [16| we analyze the attempts at extracting physical predictions for models at low 
N by computing higher orders of the 1/N expansion on finite lattices. 

In the same perspective, we discuss in Sect. [17| an alternative approach to 1/N- 
expandable spin models based on truncated Schwinger-Dyson equations. 

Sect. [18] is devoted to a review of the results that can be obtained by applying the 
methods discussed here to a wide class of 1/iV-expandable fermionic lattice models. 

Finally in Sect. [19] we briefly discuss the relevance of our results and draw our conclusions. 



2. SPIN MODELS IN d DIMENSIONS AND RENORMALIZ ABILITY 

The renormalization-group properties of two-dimensional spin models, notably asymp- 
totic freedom, are the foundations of our belief in the existence of nontrivial renormalized 
quantum field theories describing the critical behavior of these models in the coupling do- 
main lying in the neighborhood of the (trivial) critical coupling g c = (f3 c = oo). In the 
case of 1/iV-expandable spin models, the corresponding field theories are nonlinear a models 
defined on symmetric spaces. These in turn may also be thought as special limits of linear 
a models (sometimes coupled to gauge fields) enjoying the proper group symmetries. 

In order to understand properly the renormalization-group properties of these models, 
it is certainly convenient to extend the analysis by treating the physical dimension d as a 
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continuous parameter in the range 2 < d < 4. It is now possible to compare the e = d — 2 
expansion of the nonlinear models with the e' = 4 — d expansion of the linear models. These 
latter theories are known to be superrenormalizable, with an ultraviolet-stable fixed point 
at the origin and an infrared-stable fixed point at "strong" couplings when e' > 0. Since 
the infrared limit of linear models has the same relevant operator content as the ultraviolet 
limit of the nonlinear models, the latter must also be renormalizable when e < 2 around 
their nontrivial ultraviolet fixed point P4] , |T5| |. As a consequence, the critical exponents 
should take the same value when computed in the linear and nonlinear models at the same 
dimensionality, and in particular the critical exponents of the nonlinear a models should 
become trivial at d — 4. As we shall see, this phenomenon is beautifully illustrated in the 
1/N expansion. 

In the two-dimensional limit e — > however general theorems fl6l , |T7|j insure the im- 
possibility of spontaneous breaking of continuous symmetries. Therefore there is no weak- 
coupling, broken-symmetry phase, and g c = is an ultraviolet fixed point, around which 
logarithmic deviations from scaling are allowed (asymptotic freedom). Since no massless 
modes can be present, nonperturbative mass generation must occur. 

This by now standard scenario has been the starting point for most perturbative studies 
of nonlinear a models. After Polyakov's pioneering paper [IBf, perturbative ultraviolet 



renormalizability was discussed by Brezin, Zinn- Justin, and Le Guillou 



i-i 



,]19|j2g] and by 
models. The 



Bardeen, Lee, and Shrock |^TJ for O(A^) models, and by Valent |22| for CP 
extension to more general symmetric spaces was suggested by Eichenherr and Forger p3|J2~ 
and discussed by Pisarski f25j, by Duane p6|| , and by Brezin and coworkers |27[| . Quite 



naturally, one adopts dimensional regularization and evaluates the renormalization-group f3 
and 7 functions in the minimal subtraction (MS) scheme. However, a rigorous treatment 
of these models shows that generic Green's functions are plagued in two dimensions by 
severe infrared divergences. This problem was first tackled by Jevicki p8| and by Elitzur 
|29| |, who showed that two-dimensional Green's functions that are fully invariant under the 
symmetry group of the model could be computed (in low orders of perturbation theory) 
and found to be infrared-finite. This property was exploited in Refs. |3O-02]1 in the context 



of dimensional regularization, and was given a rigorous proof to all orders of perturbation 
theory by David f33^05|. Three-loop calculations for the renormalization-group functions 
were first presented by Hikami and Brezin |36[ for 0(N) models, and by Hikami |37] for 
QpAf-i moc } e } s Extensions to more general symmetric spaces was given in Refs. |38| , |39 



Anomalous dimensions were computed by Wegner and collaborators to three-loop order 
fjOfl , and later to four- loop order |41|,|42| . The four-loop order /3-function computation was 
completed in Refs. It would be beyond the purposes of the present review to give 

any detail of the abovementioned computations. A collection of results is presented in 
Appendix |A| for easy reference. 

In the context of the perturbative approach, another thoroughly-investigated issue is the 
classical thermodynamics of the models at nonzero external magnetic field, which ensures 
the absence of the infrared divergences discussed above. Renormalization and scaling behav- 

H| , who exploited 



ior were discussed in Ref. |1J] and reconsidered by Jolicoeur and Niel [|Z 
the scaling properties to devise an extrapolation method allowing for nonperturbative pre- 
dictions in the limit of vanishing magnetic field. 

When we turn to the approach based on the 1/N expansion, willing to investigate the 
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renormalizability properties of the models for 2 < d < 4, we must dramatically change our 
focus from a situation where the parameter e can be considered infinitesimal and employed 
as an ultraviolet regulator to the case where the physical dimensionality is a fixed finite 
parameter. Nothing prevents us in principle from using some of the dimensional regular- 
ization techniques, and this is one of the basic ingredients in the study of rf-dimensional 
models around criticality and in the evaluation of critical exponents in the 1/N expansion. 
An impressive series of results were obtained by Abe |47]-H9|], Brezin and Wallace [50|, and 



Ma |pl| , p2| , and more recently with improved techniques by Vasilev and coworkers |p3|-p6 
In O(N) models, the critical exponent rj is by now known to 0(1/N 3 ), while the critical 
exponent v is known to 0(1/N 2 ); 0(1/N) results are available for CP^" 1 models. Following 
the procedure indicated in Ref . |36| , these results are also the starting point for an evaluation 
of the renormalizat ion-group /3 and 7 functions (in dimensional regularization and minimal 
subtraction) at the same orders of the 1/N expansion. All the results for critical exponents 
confirm the abovementioned observations about universality between nonlinear and linear 
models and triviality in d — 4. 

We want to stress that, as far as the present evidence goes, the e and 1/N expansions 
appear to be strictly commuting when applied to the evaluation of physical quantities, such 
as the critical exponents. 

The problem of renormalizability in the framework of the 1/N expansion was studied by 
Symanzik p7||58l and by Arefeva and collaborators. In Refs. f59H63| the ultraviolet renor- 
malizability of the three-dimensional O(N) models in both the symmetric and the broken- 
symmetry phase was shown to all orders of 1/N by applying dimensional regularization. 
The ultraviolet renormalizability of CP^ -1 models when d = 2, 3 was shown in Ref. |3J by 
similar methods. However, in order to prove the existence of a renormalized critical theory 
free of infrared divergences, it was originally necessary to give up dimensional regularization 
and attack the problem from the point of view of BPHZL renormalization, which was done 
for three-dimensional O(N) models in Ref. |]B5|| . Subsequently the result was generalized to 
all 2 < d < 4 by the introduction of analytic regularization |66fl . 

A different view of renormalizability for asymptotically-free 1/iV-expandable field theo- 
ries has been put forward by Rim and Weisberger pTfl . The essential, if subtle, equivalence 
of this point of view with more standard dimensional regularization approaches has however 
been exposed in Ref. |5|]. 

A very important issue in the context of the 1/N expansion of CP N ~ 1 models is the 
relevance of classical instanton configurations, that appear to be nonperturbative in the ex- 
pansion parameter 1/N, and therefore might in principle invalidate conclusions obtained in 
a purely perturbative context. The problem was however solved by Jevicki pSf, who showed 



that, at the quantum effective- action level, instantons, instead of being stationary points, 
appear in the form of poles. One may then demonstrate that the 1/N expansion and the 
semiclassical method correspond to two alternative contour integrations of the functional in- 
tegral. Further insight on the role of instantons was obtained by David ||69|| , who discussed 
also the problem of summability of the instanton contributions, computed in Refs. |7D| , |7T| . 
The quantum statistics of CP^ -1 models was studied by Affleck [0,[f3|] also in connection 
with the topological properties (^-dependence) of the models, and later analyzed and re- 
viewed by Actor [[74[]. An extension to CP^ -1 models coupled to fermions was discussed 



in Ref. [|75||. For sake of completeness, we also mention that a different nonperturbative 
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approach to CP^ -1 models, based on the role of "torons" (classical solutions with fractional 
topological charge) has been put forward in recent years by Zhitnitsky [JFE] and found to 
agree with large- N predictions. 

A more general nonperturbative issue that may be addressed in the context of the 1/N 



expansion is the existence and the role of infrared renormalons |77| , appearing as singularities 
on the positive real axis of the Borel transform in massless ultraviolet-free theories, and 
related to the appearance of nonperturbative expectation values. These in turn are the basic 
ingredients in the operator product expansion approach advocated by Shifman, Vainshtein 
and Zakharov in order to describe large- distance effects in asymptotically free theories 



David showed that, in the context of the 1/N expansion of O(N) models, nonperturbative 
terms can be organized in an operator expansion, but they have infrared renormalons |79| ; 



these renormalons cancel against the corresponding renormalons appearing in the coefficients 
of the operator product expansion when Green's functions (involving only zero- dimension op- 
erators) are computed. According to the same author |80|j8T|] , only in well-definite instances 
(e.g. the topological charge density, and other quantities with a direct physical meaning) 
nonperturbative expectation values can be defined unambiguously. In any case, it is possible 
to show that, in each order in 1/N, the O(N) two-dimensional .S-matrix amplitudes can be 
written as series in powers of the dynamically-generated mass times a convergent perturba- 
tive series [[B^l - For the partially different point of view supported by the ITEP group, one 
should see Refs. |83|]84}| , where the issue of the operator product expansion in the context of 
the 1/N expansion of O(N) models is also discussed. 

The subject of the operator product expansion and renormalizability for critical O(N) 
models in dimension 2 < d < 4 (where nontrivial criticality exists) has been thoroughly 
investigated in recent years by Lang and Ruhl []55|-|9"0||. 

Finally we should mention that the possibility of including nonperturbative effects di- 
rectly into the perturbative expansion has been explored by Davis and Nahm, who discussed 
both 0(N) models and CP^ 1 models f92|, showing that proper normal-ordering may 



lead automatically to the inclusion of a nonperturbative mass gap in the perturbative series 
(and confinement in CP^ -1 models), and the result of this procedure commutes with the 
1/N expansion |)3]. In a related development J94j the vacuum structure of the 0(N) model 



is studied by a variational technique and agreement with conventional large-iV results is 
found. 



3. 1 /TV-EXPANDABLE TWO-DIMENSIONAL SPIN MODELS 

In order to achieve some generality, we shall investigate the properties of a two-parameter 
class of 1/iV-expandable spin models, described by the continuum action 

S = N J d 2 x {PydpZdpZ + P g D^D^} , (3.1) 

where z is an iV-component complex field subject to the constraint 

zz = 1 (3.2a) 
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and a covariant derivative = 8^ + iA^ has been defined in terms of the composite gauge 
fields 

= \i {zd^z - d^z z) = izd^z. (3.2b) 

This action was introduced first by Samuel and it is an interpolating action between pure 
CP^ -1 models (/? v = 0) and U(iV) vector models (/3 g = 0), which in turn are nothing but 
0(2iV) vector models. 

We notice that pure CP^ -1 models enjoy a U(l) gauge invariance related to the local 
transformations 

z{x) -> e iX{x) z{x), (3.3a) 
z{x) -> e- iX ^z(x), (3.3b) 
A^(x) -> A M (x) - 9 M A(x). (3.3c) 

This invariance will play an important role in determining the structure of the effective 
vertices. 

We can introduce the (rescaled) weak coupling parameter / by the following change of 
variables: 

^ = /3v + /3 g , (3.4a) 

*/ = (r • (3 - 4b) 

The coupling constant x, as defined by Eqs. Q3.4Q , enjoys the property of non-renormali- 
zation, i.e. the renormalization-group trajectories in the (/3 S , j3 v ) plane (in the continuum 
version of the model) are just the curves of constant x and correspond to physically different 
theories. This property will emerge rather clearly from the discussion of the 1/N expansion. 
The renormalization-group trajectories are plotted in Fig. [T] for several values of x. 

The 1/N expansion is achieved as usual by implementing the constraints ( |3.2|) by 
Lagrange-multiplier fields a and 9' . The manipulations are quite standard and we obtain 

S = t^t J d 2 x ^dpzduz + ^ d » z f + ia{zz-l)+ ft ~ ^ d ^ z )j 

= Yfj ^{t^^^ + y^K^ + ^^I' + ^^-i)}- (3-5) 

We can now perform the Gaussian integration over the z fields and obtain the effective 
action 

S cS = iVTrln + ~\~^>cf + + ^(-<a), (3.6) 

where now = + i9 ' . Finally, rescaling the multiplier field 9' to 9^ = 9' /{I + x/) and 
introducing the vacuum expectation value of the a field, a(x) = (a) + a q (x), (a) ■ 
we obtain the following form of the effective action: 



-im, 
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S cS = NTr In - i 0„} + m 2 + ia q } + — {-m 2 Q - ia q + (1 + x/)^} . (3.7) 

In the large- N limit, the value of wig is determined, as a function of / only, by the saddle- 
point condition (gap equation) 

-i= r d2p 1 f38i 

2f J (2vr)2p2 + m 2- I • ) 

Eq. ( |3.8|) is in need of ultraviolet regularization; we shall come to this point in Sect. [5]. 

By taking the second functional derivative of the effective action around the saddle point 
we may now obtain the propagators of the quantum fluctuations associated with the fields 
a q and 6^; both are 0(1/ N) quantities that can be expressed by the functions 

A -i n= f^_^ i i 1 e + i r . q ^ 

^ {P> J (2tt) 2 q 2 + m 2 Q (p + q) 2 + m 2 2tt P ^ £ ~ 1 ' 1 ' 

1\ . f d 2 q (p fl + 2q ll )(p u + 2q l/ ) 



(2vr) 2 [g 2 + mg][(p + g) 2 + mg] 



x5 M , + — ( £ln|— r - 2 ) ( V - ™£ ) , (3.9b) 



where £ = a/1 + 4mg/p 2 . The effective propagators in d dimensions are presented for 
reference in Appendix [Bj 

It is now crucial to observe that all the high-order effective vertices resulting from 
Eq. (|3.7| ) by taking higher functional derivatives are completely unaffected by the value 
of x. As a consequence, the vertices share the gauge properties enjoyed by the CP^^ 1 
model, and in particular transversality. This fact in turn implies the possibility of replacing 
the propagator (|3.9|b) 

A (e)fiu (p) = A (e) (p) (V - P ~f) + ~ ^ , (3-lOa) 

where 

A ( 9 )(P) = ^ (^ ln " 2 + 27rx ) ' ( 3 ' 10b ) 

with its transverse part, when computing expectations of gauge-invariant operators. 
Eq. ( |3.10| ) shows that x is a "physical" parameter, related to the ratio of the mass of 
the Lagrangian field z (m in the large-iV limit) to the mass m# of the propagating field 
6^. The exact relationship between the two masses is expressed in the large- N limit by the 
equation 



' 4?7z 2 / Atti 2 ! 

— ^-larccotW — £ - 1 = 1 - ttx, 0<x<-. (3.11) 
me V m e 71 

This is the basic physical reason why x is not subject to renormalization. mg(x) is plotted 
in Fig. |. 

Let us now briefly describe the classes of different physical theories parametrized by x. 
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1. When x = we get m« = 0: 9^ becomes a dynamical gauge field giving rise to a linear 
confining potential between z and z, and the physical states are the bound states that 
are singlets under gauge transformations. This is the well-known physical picture of 
the CP^ -1 models |95|-07l. 



2. For very small values of x the above picture is substantially unchanged: in the absence 
of a gauge symmetry the z fields are not automatically confined by Elitzur's theorem; 
however, their mass is so much higher than that of their bound states to make them 
effectively disappear from the physical spectrum. 

3. With growing x an inversion occurs and the z fields become the fundamental states 
of the model, while the mass of the bound states becomes bigger and bigger, when 
measured in units of mo. 

4. At x = 1/tt we meet a threshold: mg = 2mo and the Yukawa potential that was 
the remnant of the linear confining potential completely disappears. This is a quite 
interesting model: it is easy to get convinced that the corresponding action is nothing 
but the effective action resulting from the functional integration over a set of iV- 
component massless fermion fields minimally coupled to the "gauge" field (cf. e.g. 
Refs. 0,0). This model has no quantum anomaly, and its factorized S- matrix is 
therefore known: the physical states are in the fundamental representation of U(iV) 
and the bound states have disappeared, as expected. 

5. When x > 1/tt, there are no bound states, and the models interpolate smoothly from 
U(iV) to 0(2N) symmetry. In particular, for integer values of n = 7rx the models 



describe the minimal gauge-invariant coupling of n "flavors" of massless fermions p9 
This picture gives further support to the notion that x is a physical parameter not 
subject to renormalization. 

Finally, when x — > oo the effective field 9^ completely decouples and we are left, as 
expected, with the well-known 0(2iV) nonlinear sigma model, possessing a factorized 
S-matrix for the fundamental 2N real fields and showing absence of bound states. 



4. REVIEW OF EXACT RESULTS 

In Sect. [| we mentioned that, for special values of the parameter x, a number of ex- 
act results are available, especially concerning exact factorized S-matrices and bound-state 
spectra. We must however keep in mind that these results have been obtained in a rather 
indirect way, by applying such methods as analytic S'-matrix theory or Bethe Ansatz. As 
a consequence, the 1/N expansion offers the possibility of verifying the applicability of the 
abovementioned methods to the models at hand and therefore the correctness of the physical 
interpretation. For future reference and comparison, we would like to present here a short 
review of these exact results. 

We start from the observation that the key ingredient for the possibility of constructing 
an exact S'-matrix is the assumption of factorization of multiparticle amplitudes into two- 
particle amplitudes, i.e. absence of particle production ||100|| . This property is in turn related 
to the existence of higher-order conservation laws. Consider the standard Noether current 

associated with the global symmetry, and notice that the nonlocal charge [|101[ 
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g(m) = i y dxidx 2 e(xi-x 2 ) b'o(Mi) j (*,^2)] + y dxji(t,x), (4.1) 

where e(x) is the sign function, is classically conserved owing to the equations of motion. 
At the quantum level one may show that this conservation law is in general spoiled by 
anomalies generated by the renormalization process. These anomalies are nonperturbative, 
but their coefficients can be calculated in perturbation theory. A detailed analysis shows 
that in 0(2iV) models the anomaly is actually absent |[L02|| , while in the pure CP^ -1 case 
one finds If03| , g0l 



-^— = - / e^ l/ (t > ar)dar. (4.2) 
di 7r J 

However, when we consider the inclusion of minimally-coupled massless fermions, we realize 
that the nonlocal charge is classically conserved only if we include a contribution from the 
fermionic axial vector current. The axial current in turn is known to possess a quantum 
anomaly, opposite in sign to the r.h.s. of Eq. (|4.2j ). As a consequence, the modified nonlocal 
charge has no net quantum anomaly and the corresponding model turns out to possess a 



factorized S'-matrix ||104| , |105|| . In our language, these results imply the possibility of finding 
explicit S-matrices in the cases x = oo and x = 1 /tt respectively. 

Without belaboring on the techniques used in order to solve the factorization equations, 
we only recall that the Fock space is decomposed irreducibly into subspaces labeled by a 
definite particle number n. States are identified by sets of particle momenta {P r }, and 
factorization is expressed by 

(P' (out) | P (in)) = (p' (in) | J] S{P' r , P.) | P (in)), (4.3) 

l<r<s<n 

where S(P^, P s ) is the two-particle S'-matrix. S is best expressed in terms of rapidity vari- 
ables 9 r = arctanh(P rj i/P ri o). Imposing 0(2iV) symmetry ||106|| one may assume 



(O'iJ;02J ( out ) I u k;9 2 ,l (in)) 
= 5(9, - 9[) 5(9 2 - 9' 2 )[S ij 5 kl 8,(9) + 5*6^ S 2 (9) + 8 a 8 jk S 3 ( 
+ 5(9, - 9' 2 ) 5(9 2 - 9' 1 )[5 ij 5 kl S,(-9) + S u S jk S 2 (-9) + 5 ik 5 jt S 3 (-9)] (4.4) 

where 9 = 61 — 6 2 . Solving the constraints one finds 

S 3 (0) = - S 2 (6), (4.5a) 

s ^ = - { N-mn-e) s ^ (4 - 5b) 

and, assuming minimality in the number of poles and zeros in the physical sheet, 

S 2 (9) = R(9)R(itt -9), (4.6) 

where 
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r( , 1 *Wl-» 



; 2 2(iV- 1) 2vr J V 27r / 



The result shows no bound state poles ||107||, and was checked in a 1/N expansion up to 



second order [|108|| . 

Qpjv-i moc [ e } s w ith minimally-coupled fermions in turn correspond to SU(iV) symmetry, 



and the form of the factorized S- matrix is [|9£, 109| , 110 



(9[, i; 9' 2 ,j (out) | 9 U k; 6 2 , 1 (in)) = 5(9[ - 9,) 5{9' 2 - 9 2 )[5 ik 8 fl u x {e) + 5 jk 8 u u 2 { 

- 5(9[ - 9 2 ) 5(9' 2 - 9i)[6 jk 6u u x {9) + 5 ik 5 jlU2 (9)}, (4.8a) 

(9' 1 ,i; 9' 2 , j (out) | 9 U k] 9 2 , 1 (in)) = 5(9[ - 9 X ) 5(9' 2 - 9 2 ) [5 lk 5 3l t x (9) + <f w <f - t 2 (0)] 

+ - 2 ) - 9 1 )[5 ik 6 jl r x {9) + 5 kl 5 l3 r 2 {9)l (4.8b) 

where the bar indicates antiparticles. The constraints imply 

2?7T 

u 2 (9) = -— Ul (9), (4.9a) 

n{9) = r 2 {9) = 0, (4.9b) 

t 2 W--^|^yt lW . (4.9c) 

Moreover, the crossing symmetry requires 

ti(0) =ui(z7r-0), (4.9d) 

* 2 (0) =u 2 {iir-9). (4.9e) 

Finally, from minimality one obtains 

1 

'■ W= Jl «wi 1 tfV (410) 





) r 




1 

-iv + 


Z0 


2vry 


s- 


2^ 
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Z0 


2tt j 






2^ 



2 

There are no bound state poles: the bosons interact repulsively and the fermions are screened 
by a "secret" long-range force, while the gauge field loses the zero mass pole. This result 
was again checked in the 1/N expansion. 

We must notice that these integrable models can sometimes be solved by the Bethe 
Ansatz and quantum inverse-scattering methods, always reproducing the abovementioned 
results. 

For completeness we mention that the factorized 5-matrix approach allows in some spe- 
cial cases the determination of exact form factors. In particular, in 0(2iV) models it is 
possible to evaluate the matrix elements of the Noether current between the two-particle 
state and the vacuum. The result was checked in the 1/N expansion to 0(1/N) by the use 
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of the (explicitly known) spectral representation of the propagator A^(p 2 ) ||111|| . One can 
write 



sinh ib 

.- \ i — 



A( a )(p 2 ) = Anm\ 



sinh ^ip 



Ami 



(4.11) 



Therefore 



A( a )(p 2 ) = Anml 



1 + 



V 

2^ 



d ^ 2 P(a)(/i 2 ) 
!m 2 fi 2 (p 2 + ^ 2 ) 



(4.12) 



2vr 



sinh0 

</> 2 + 7T 2 ' 



cosh 



2 1 



p 



Ami 



Likewise, when k— 1/tt we have 



A*)(p 2 



2tt 



tanh ^ 



p 2 +p 2 



d/i 2 



(4.13) 



pw(p 2 



2vr 



coth i( 



+ 7T 2 



A last very important exact result that can be obtained from the analysis of the S- 
matrices for integrable models is the analytic determination of the so-called mass - A- 
parameter ratio, where the A parameter is defined in standard perturbation theory (MS 
scheme) in terms of the subtraction scale fi and the universal coefficients bo and h\ of the 
renormalization-group (3 function: 



^MS = P eX P 



dg' 



-bi/bl 



exp 



1 

bog 



(4-14) 



The seminal result was obtained by Hasenfratz and collaborators in the case of 0(2N) 
models [|112jJ113|| : 



m 



8\ 2(AT-1) 

e 



1 



r i 



A 



MS • 



2(N- 1) 

The analysis can be repeated in the x = 1 /n model; the final result is 



m 



j_ 

2\ n 

e 



ni+ N 



| \ ^MS • 



(4.15) 



(4.16) 



As we shall see in the next Section, these results can also be explicitly verified in the context 
of the 1/N expansion. 



15 



5. THE 1/JV EXPANSION IN THE CONTINUUM: REGULARIZATION AND 

RENORMALIZATION 



The Feynman rules for the 1/N expansion of the model can be easily derived by intro- 
ducing external currents coupled to the z fields before performing the functional integration. 
They are summarized in Fig. |3]. The 1/N expansion is an expansion in the loops of the effec- 
tive fields a q and 9^. In the graphical representation, a closed loop of the z fields stands for 
an effective vertex of the expansion; effective vertices can be obtained by taking functional 
derivatives of the effective action and carry a factor of N. 

The effective vertices actually amount to one-loop integrals over the fundamental field 
propagators. In two dimensions, all the effective vertices may be in principle computed 
analytically, but the computation may become very cumbersome in the case of exceptional 
configurations of momenta, which are often those relevant to the actual computations one 
would like to perform. For a discussion of this technical problem, cf. Ref. [114|. 



While no regularization is needed in the evaluation of the effective vertices, it becomes 
unavoidable when one wants to compute the Green's functions of the physical fields. In our 
choice of regularization we were not guided by the usual requirements of Poincare invari- 
ance, gauge invariance, and consistency up to all perturbative orders that made dimensional 
regularization a favorite tool of quantum field theory. Having in mind our final purpose 
of performing explicit lattice computations, we rather focused on the requests of computa- 
tional ease and transparency in the regularization mechanism, which is often obscured in 
dimensional regularization by the interplay of ultraviolet and infrared singularities. 

Relaxing the consistency request down to one-loop consistency, we found that the sim- 
plest and most transparent scheme was a kind of sharp-momentum (SM) cutoff. The regu- 
larization procedure (roughly formulated for O(N) models in Ref. ||115fl and discussed more 
precisely in Ref. |||) starts from the observation that, since the A propagators are finite, the 
0(1/N) (one- loop) Feynman integrals appearing in the computation of the n-point Green's 
functions are integrals of regular functions. They can therefore be regularized by subtracting 
explicitly the highest powers of the integration variable appearing in the Taylor expansion 
of the integrand. The lower limit for the integration of the subtraction terms is arbitrary, 
and we are therefore introducing a dependence on the cutoff value M 2 . When the integrals 
are finite, this dependence disappears when taking the limit M 2 — > oo. 

Let us however consider the regularized gap equation 

1 f d 2 p 1 f°° dp 2 1 1 , M 2 . . 



2f J (2tt) 2 p 2 + m\ J M 2 Air p 2 Air m 

Eq. (|5.1|) allows us to eliminate the dependence on M 2 of any superficially divergent di- 
agram in favor of an explicit dependence on the coupling constant, which in turn will be 
reabsorbed in the renormalization-group-invariant definition of the physical mass and in the 
wavefunction renormalization. 

A first obvious analogy with the lattice formulation lies in the fact that we are working 
with a "bare" coupling constant which can be varied together with the cutoff while keeping 
all physical quantities constant. As we shall show in Sect. |9|, the analogy can be made much 
more stringent by finding the relationship between the SM cutoff and the lattice cutoff, 
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which will also lead to a natural regularization of the infrared singularities appearing in 
massless lattice integrals. 

It is important to observe that our regularization procedure differs by 0(1/M 2 ) terms 
from a naive sharp-momentum cutoff (i.e. simply setting the upper integration limit to M 2 ). 
This difference is irrelevant in the continuum, but it will become crucial when discussing the 
relationship with the lattice scheme. 

The connection between SM regularization and dimensional regularization (and renor- 
malization) has been explored in some detail in the asymptotic regime of large Euclidean 
momenta, where perturbation theory holds We shall not belabor on this point in the 
present paper. 



A. The free energy 



SM regularization leads to a simple parametrization of the perturbative tails that con- 
tribute to the nonscaling part of physical quantities. We can therefore evaluate the scaling 
contributions to the free energy to 0(1/N) in our models, thus also offering a first explicit 
example of our computational techniques. The free energy is the sum of the connected vac- 
uum diagrams of the effective theory. In lowest orders only the trivial Gaussian integrations 
over the z, a and 6^ fields do contribute. Subtracting the perturbative tail and keeping 
only the scaling part, according to the rules of SM regularization, leads to the following 
expression for the first two nontrivial contributions: 



F = A Tr In 



where 



p 2 + ml 
p2 



A {a)ip) 

A*)(p) 



A 



mi H — Tr In 
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(P) 



+ 



2-rcp 2 



p^oo \n(p 2 /ml) 



A® 

(a) 



(p), 



2tt 



p^oo 



A {0) 



(P)- 



ln(p 2 / m o) - 2 + 27rx 
According to our rules, the regularized expression of the free energy is therefore 



4tt 2 



d 2 p 



(27T) 
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£-1 

1 2ttx 
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dp 2 2ml 



3 — 27TX 



M 2 An p 2 \n(p 2 /ml) + 2nx — 2 



In 



+ 



M 2 
p 2 

In — k + 2ixk 
1 

A 



ln(p 2 /ml) 
-2 



47r 4ir 



M 2 I M" 

In In — + (3-2ttx) In I In — + 2tt 



M 2 



2 + c F (x) 



(5.2) 

(5.3a) 
(5.3b) 



O ^ 



1 

A 



(5.4) 



where cp(x) is a numerical constant which can be computed to all desired numerical accu- 
racy; it is plotted in Fig. [|. Notable special cases are 
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c F (0) = 1.18887122, 
1 



c F 



7T 



2 7 ; 



E- 



(5.5a) 
(5.5b) 



The issue of the evaluation of dimensionless SM-regulated one-loop continuum integrals is 
discussed in Appendix |C| 

Cp(x) may also be evaluated in the context of a 1/x expansion of Eq. ( |5.4| ). The result 

is 



c F {x) = (27rx-3)ln2vrx + 7E + ln4-3 + — - + 0[ ) . 



The large- x limit of Eq. (|5.4j) is [116 



2 9 

mn mt 



F = N - H - 

A.TT 4:71 



2ir x 



M 2 , M 2 

In In — 5- + 7e — In - — - — 1 



4m 2 , 



Substituting Eq. ( pM|) in (5.4) we obtain 



F = N — - H 

47T 47T 



2tt 



./ 



2ty 



In — + (3 - 2vrx) In — + 2ttx - 2 + c F (x 



./ 



+ 



N 



(5.6) 



(5.7) 



(5.8) 



Now, exploiting the renormalization-group invariance of the scaling part of the free energy, 
we are ready to obtain the renormalization-group (3 function: 



/?(/) 



d\nF 



7T 



Iff, 3-27TX 
1 + TTf- 1 + 



N 2ir \ l + /(x- 1/tt) 



f-L 

\N 2 



(5.9) 



It is easy to check that the known universal coefficients are correctly reproduced. In partic- 
ular, when x = 1/tt 



P(f) = ~ 



f 2 



IT 



1 + 



11 

N 7T 



+ 



and when x — > oo 



N 2n 



1 

iV 2 



(5.10) 



(5.11) 



as expected for 0(2iV) nonlinear sigma models. 



B. The two-point function: regularization 

We now focus our attention on the properties of the invariant two-point correlation 
function of the fields in the fundamental representation. Let us define 

G(P) = i / d 2 xe^ (z(x)z(0)) = -^— 2 - i — Uf27i(p) — J-^ +o(±) ■ 
2 j J p 2 + tuq N p 2 + rriQ p 2 + m,Q \N 2 J 

(5.12) 
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The 0(1/ N) contributions to the two-point function are drawn in Fig. [5[ We obtain 



Slip) 
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+ ml 



(5.13) 



By straightforward manipulations (essentially replacing zero-momentum insertions of the a q 
field with derivatives with respect to m^), Eq. ( |5.13| ) can be recast into the form 
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(5.14) 



Furthermore, explicit knowledge of the propagators allows us to make use of the identities 
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(5.15a) 
(5.15b) 



As a consequence, we obtain the representation 
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(5.16) 



The last term in Eq. (|5.16|) reflects the dependence on the longitudinal degrees of freedom 
of the field. Its contribution can be computed in closed form (after regularization) and it 
amounts to 
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— (p + m ) 
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1 , M 2 
— m — 7T 
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1 ln P 2 + m 2 
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2tt 
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(5.17) 



Eq. (|5.17 ) is singular in the x — > limit, when the model becomes gauge- invariant and the 
contribution of the longitudinal degrees of freedom becomes gauge- dependent. 

The regularized version of Eq. ( |5.16|) is obtained by applying the SM scheme prescrip- 
tions: 
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Eq. (|5.18|) implies the possibility of parametrizing El eg (p) in the form 

r m 2 / m 2 \ 
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(5.18) 



(5.19) 



where Uf n (p) is a regular, M-independent function of p 2 (and x) subject to the normalization 
condition 



Z™(im Q ) = 0, 

and c m (x) is a numerically computable constant. 



(5.20) 



C. Mass and wavefunction renormalization 



The interpretation of Eqs. (|5.19| ) and (|5.20 ) becomes straightforward in the context of 
renormalization, when we identify G(p) with the bare two-point function and write 



1 



p 2 + m 2 + 1 Z? g {p) p* + m 2 + 1 £f*(p) 



(5.21) 



Eq. ( |5.20|) is then the on-shell renormalization condition and 
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(5.22a) 
(5.22b) 



In particular Eq. ( p.2C| ) insures that the mass gap m? is identified as the pole of the two-point 
function and its 0(1/N) contribution m\ satisfies the condition 



m\ = Sl eg (imo). 

Applying Eq. ( |5.23| ) directly to Eq. (" |5 . 1 8|) we obtain the representation 
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(5.24) 
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allowing a direct determination of c m (x), which is plotted in Fig. |4j. 

In the x — > limit, c m (x) is an infrared-divergent quantity, which shows that in the 
CP^" 1 models no single-particle mass for the fundamental fields can be consistently defined 
This is a reflection of the gauge properties of the model and is further evidence 
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for confinement. More generally, it is pleasant to notice that m\ does not depend on the 
longitudinal degrees of freedom of the vector field, as a consequence of gauge-invariance 
of the couplings and the on-shell condition. This allows a direct physical interpretation 
of Eq. ( |5.24|) also in those instances (e.g. the x = 1/n model) where we are dealing with 
a gauge-fixed version of a theory enjoying gauge symmetry properties. Actually we can 
compute exactly 



.(1/0 = 2 7E 



(5.25) 



Eq. ( |5.25| ) can be shown to be consistent with the 1/N expansion of the exact result ( |4.16 ) 
The 1/x expansion of Eq. ( |5.24| ) leads to 
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The large- x limit of Eq. (ggg ) is jTT6HTT8 
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consistently with the exact result ( |4.15|) . In the same limit one may also show that 
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(5.28) 



and the on-shell renormalization condition can be imposed on Z. 

It is pleasant to notice that Eq. (|5.22| a) allows an independent determination of the 
renormalization-group f3 function and the result is completely consistent with Eq. ( |5.9|) . As 
a consequence the adimensional ratio 



(5.29) 



is universal and scheme-independent. It is interesting to notice that one obtains from 
Eqs. ( |5.6| ) and ( 5.26|) the relationship 



c F {x) - c m (x) 



(5.30) 



An alternative renormalization-group-invariant definition of the correlation length can 
be defined starting from the second moment of the two-point correlation function 
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In momentum space this definition leads to the relationship 
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Substituting Eq. (|5.19|) in Eq. (|5.32|) and comparing with Eq. (|5.22| a), we can easily show 



that 



m 2 R = m 2 + 5m 
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(5.34) 



Keeping also in mind Eq. (|5.20| ) , we come to the conclusion that 5m R is amenable to a (typi- 
cally small) calculable constant, which can be interpreted as a universal scheme-independent 
adimensional ratio. The 0(1/N) contribution to 5m R is plotted as a function of x in Fig. |6|. 

Eq. ( |5.32|) becomes singular in the x — ► limit, when only gauge-invariant correlations 
can be sensibly defined. The strategy for such computations is discussed in Ref. | |114| | and 
is not especially relevant to our present analysis. We only mention that gauge invariance is 
obtained at the price of introducing "strings" connecting the z fields, and defining 



Gc{x,y) 



v ( ; 



z{x)exp\i / dtJJt)\z(0) 



(5.35) 



String renormalization is required. It is then possible to define (x 2 ) c in analogy with 
Eq. (£TT]). 

We stress that the gauge dependence of the definition Q5.31D makes the observable m R 
dependent on the longitudinal degrees of freedom of the vector field. Therefore it has no 
direct physical interpretation in those special cases when we want to recover the gauge- 
invariant properties of a gauge-fixed model. Similar considerations hold for the so-called 
magnetic susceptibility 



X = d 2 x(z(x)z(0))=2fG(0) 



ml N rriQ 



(5.36) 



The value of E[ ee (0) can be computed analytically for special values of x; ignoring the 
contributions of the longitudinal degrees of freedom we have (cf. Ref. | 116| ) 
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(5.37a) 
(5.37b) 



where 
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ci = ll /(j)r(j) ^ 
1 ni)r(l) 



0.4861007. 



(5.38) 



We can adopt a wavefunction renormalization condition slightly different from 
Eq. ( |5.22| b), defining the renormalization constant Z by 



vm| 
2/ 



1 



l--Z[(0) + O[ — ). 



iv 



(5.39) 



It is easy to check that the ratio of the two definitions is a /3-independent constant: 



f-L 



(5.40) 



Finally it is possible to compute the renormalization-group function 7(/), after its defi- 
nition 



-r(f) = -P(f)^H2fz{f)]. 



(5.41) 



Ignoring the contributions from the longitudinal degrees of freedom of the vector field, we 
obtain 



1(f) 



f 



TV 



1 - 



1 

2N 



1 



LL 

N 7T 



1 



1 



1 + /(x-1/tt) 



O 



(± 



(5.42) 



and, specifically for 0(2N) models. 

./ 



7(/) 



7T 



2N 



1 + 



LL 

N 2rr 



(5.43) 



Here and in the following, when checking agreement with the expressions presented in Ap- 
pendix [A], one must perform an appropriate change of variables, whose form may be extracted 
by comparing SM and MS j3 functions. 



D. Correlations of composite operators 

Another very important class of correlation functions is obtained by considering the 
Green's function of the (gauge-invariant) composite operators 

Pij{x) = Zi(x) Zj (x). (5.44) 

We have shown in Ref. ||114|| that, for CP^ -1 models (x = 0) 

B(x — y) ( 1 \ 

G ijtk i(x -y) = (Pij(x)P kl (y)) coiin = N ^ N + ^ Uafyfc - ^ $iM J » ( 5 - 45 ) 

where 
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B(x - y) 



NjziWzjWziMzjjy))-! _ Af 2 
N - 1 N 



(5.46) 



and in turn A 1 = NA,, and A 1 1 is expressed to 0(1/N) by the sum of Feynman diagrams 



drawn in Fig. [?]. More generally, A^ 1 is the sum of all the diagrams such that the external 
a legs emerge from the same effective vertex (generalized tadpole contributions). 

This analysis applies almost literally to the more general case presented here, and we can 
quote the final result which is just a very slight generalization of Eq. (8.20) in Ref. [ 114|| : 
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A {a) {k) 



v} a \ P , k) + vr(p, -k) + vrtp, k) + 



r (a) 



(&), 



d 



k 2 + Ami d m o 



A {e) (k)Uk 2 + Am 2 ) Vt\p,k) + vr(p,-k) - 4 [V 3 (p, k) + V 3 (p, —k)] 



(a)/ 



+ (k 2 + Am 2 + 2 P 2 )v} b \p, k) + ^^m(p) 



k 2 + Ami d m o 



(5.47) 



The effective vertices entering Eq. (|5.47|) are drawn in Fig. g. The formal definitions of V3, 



v} a \ and V^ 0> , together with their explicit expressions in terms of elementary functions, 
can be found in Appendix [D| We stress that algebraic manipulations, based on the gauge 
invariance of the effective vertices, lead to the possibility of replacing all explicit dependence 
on the mixed vector-scalar vertices appearing in Fig. [7] with purely scalar vertices. 

Regularization and renormalization are straightforward along the lines presented in 
Ref. [114 1 . One may analyze the large k behavior of the effective vertices and find that 
the ultraviolet divergence of Eq. (|5.47|) is regulated by the SM counterterm 
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dm 2 ' 



(5.48) 



Eq. ( |5.48| ) shows that renormalized Green's functions are obtained by mass and wavefunc- 
tion renormalization, and in particular mass renormalization is once more consistent with 
Eq. (|5.4|) and Eq. ( |5.22| a), while wavefunction renormalization is obtained by defining 



2.n y -21n( 7 



2ttx- 2 



(5.49) 



We can extract from Eq. (|5.49|) the anomalous dimension of the composite field in the 
SM regularization scheme: 



Mf) 



- Pif) ^ln[Af 2 Z P ] = 2 -l 



N 2tt 



1 + 



l + /(x-l/7r)yj 



o 



1 

N 2 



(5.50) 



We notice that, in contrast with Eq. (|5.22[ b) and consistently with the gauge properties of 
Pij, Zp is independent of the longitudinal degrees of freedom of M . 
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A magnetic susceptibility and a second moment of the correlation function can be defined 
for the field P^-. These quantities can be shown to satisfy all the renormalizat ion-group 
requirements. In particular 



Jd 2 x\x 2 (tr{P(x)P(0)}> 



(5.51) 



' p f d 2 x(tr{P(x)P(0)}) 
can be used as an alternative definition of the correlation length. In the large- N limit, 



'o 



In analogy with Eq. ( |5.39| ), we can obtain a computationally convenient definition of 
renormalization constant Zp by the prescription 



2tt xp 



2tt 
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dp 2 



-i 



p 2 =0 



where 



and we fixed the normalization by noticing that 

-i -l 



Xp = / d 2 x(tr{P(a;)P(0)}) 
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(5.54) 
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dp 



(5.55) 



p 2 =0 



In order to regularize Eq. ( |5.53| ), we apply Eq. ( |5.48| ) and recognize that the structure of the 
counterterms of A 2 sAj is simply 
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+ 2(3 - 2ttx) 



(5.56) 



_d_ 

dp 2 



(5.57) 



p 2 =0 



we immediately check that Eq. ( |5.53| ) is consistent with the parametrization ( |5.49| ); cz{>c) 
can be computed numerically and is plotted in Fig. |[ It is worth noticing that c^(x) is 
finite for k —>■ 0. 
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E. Wilson loops and static potential 



Having exhausted the discussion of the correlations of fundamental fields that may be 
relevant to a 0(1/ N) analysis, we would like to consider also the properties of (gauge- 
invariant) correlations of the vector field 9^. These correlations are most naturally expressed 
in terms of expectation values of the Wilson loops. The vector field being Abelian, no path 
ordering is required and the general definition, for arbitrary loops C, is 



L(C) = ^exp ji jf dt M M (t) 



1 " <b dt u & dt[ 



2N 



d 2 k 
(2tt) 2 



(5.58) 



We defined L(C) having in mind the interpretation of our Lagrangian as an effective theory 
for an underlying gauge-invariant model where the gauge field is 9^, not the original vector 
field 9' 

For our purposes we shall only consider long rectangular loops of size RxT, and consider 
the limit T — > oo. The quantity 



V{R) = - lim — In L(R, T) 

T — »oo ± 



(5.59) 



can be interpreted as the interaction potential generated by vector fields between two static 
sources. This quantity is relevant to the discussion of the nonrelativistic bound-state spec- 
trum, a reasonable approximation in the large- N limit. One can easily show that V(R) does 
not depend on the longitudinal components of 9^, and within the 1/N expansion one may 
generalize the result of Ref. [114] to 
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NV{R) ^ -2 / cos kR A {e) (k) — + 2 
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dk 

J 2^ 



(5.60) 



where A^(k) is defined in Eq. (p.3|b); it has been introduced in the context of SM regu- 
larization of the loop ultraviolet singularities giving rise to the so-called "perimeter term" . 
Rotating the integration contour in the complex k plane to k = ix + e we finally obtain the 
following representation of the static potential, for < x < 1/tt: 
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27TX 
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7T 2 £ 



m 



(5.61) 



where £' = y 1 — Arn^/x 2 and we have defined the mass nig according to Eq. (|3.11|) . In 
Fig. [10] we have drawn the function NV(R)/mo for a few different values of x. 

Eq. ( |5.61|) , considered as a function of x, interpolates between the linear confining po- 
tential at x = and the limiting case x = 1/tt, when lim^oo V(R) = and the bound state 
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spectrum finally disappears, while me = 2m . The qualitative discussion of the properties 
of the models as a function of x, presented in Sect. [| is essentially based upon an analysis 
of this static potential, plus the information coming from integrability at x = 1/ir. 

A more quantitative discussion of the bound state spectrum can be obtained, for 
sufficiently large N, by considering the semiclassical approximation, i.e. by solving the 
Schrodinger equation in the asymptotic potential Vas(R). Ks(-R) is obtained by remov- 
ing from Eq. ( |5.61| ) all the contributions not affecting leading order predictions, and it can 



be represented in the form 



VWfl,*) = ^4^)(l-e-™.«), (5.62) 



where 



x 



1 - -x 2 

A(x) = -- 4 (5.63) 

6 \x 2 - 1 + - 1 arccot y^A/x 2 - 1 

enjoying the property 

\imxA(x) = 1. (5.64) 

x~ »0 

The analysis of the resulting Schrodinger equation is presented in some detail in Appendix [E]. 
In the case x = (CP^ -1 models), the analysis of the Schrodinger equation in the linear 
potential was presented in Refs. |97|J119|| and in more detail in Ref. [|93|]. 



F. Topological charge and susceptibility 

Before concluding this Section, we must mention that, in the special case x = 0, the 
1/N expansion can also be applied to the problem of the so-called topological susceptibility 



X t = j d 2 x(q(x)q(0)}, (5.65) 

where 

1 



q(x) = — DpZ D v z = — e^dffiy (5.66) 



is the topological charge density. One can easily show that 

1 



Xt = hm P A ($) (p), (5.67) 



v 



where A^(p) is the full propagator of the quantum field 9^. Xt is trivially zero when x ^ 0. 
In the large iV limit, Eq. ( p.lOj b) implies the simple relationship 
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The computation of the 1/N 2 corrections to the inverse vector field propagator was per- 
formed in Ref. [ 120|| . They can be obtained from the diagrams of Fig. 11. The vector and 
mixed scalar-vector vertices can be replaced by combination of scalar vertices; the result is 
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(5.69) 



where 



W!(k,p) = -(p 2 + ±ml)[Vt\k,v) + V} a \k, -p)) 

- (p 2 + 2k 2 + 4m 2 )v} b) (k,p) + A[V 3 (k,p) + V 3 (k, -p)], 
W 2 (k,p) = (p 2 + Am 2 )(k 2 + Am 2 )[v} a) (k,p) + v} a) (k, -p)] 
+ (p 2 + 2k 2 + Am 2 )(k 2 + 2p 2 + Am 2 )v} b) (k,p) 
- A(p 2 + k 2 + Am 2 )[V 3 (k,p) + V 3 {k, -p)} 
k 2 p 2 



(5.70a) 



(k-p) 



:[A {a) (k + p)Zi(k,p) + A {a) (k-p)Z 2 _(k, : 



and in turn 



Z±(k,p) = (p 2 + k 2 ± Am 2 )V 3 (Tk,p) - A^(k) - A^(p). 



(5.70b) 



(5.71) 



A property of W 2 (k,p) relevant to the computation of higher-order corrections to the slope 
of the linear static potential is 



r W 2 (k,p) 
lim /2 2 

p^o k z p z 

k-*0 
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457r m! 



(5.72) 



o 



Ultraviolet regularization in the SM scheme is straightforward and can be proven to be 
consistent with the expected renormalizability properties of the model: only mass and z- 
field wavefunction renormalization are required. 

The 1/iV-expandable, dimensionless ratio R = Xt (% 2 )p is computed; it takes the scheme- 
independent value 



R 



2nN 



N 



N 2 



(5.73) 



where c R ^ -0.380088. 



6. LATTICE FORMULATION OF TWO-DIMENSIONAL SPIN MODELS 

It is well known that infinitely many lattice actions share the same naive continuum limit. 
However, from the point of view of numerical simulations, the choice of a lattice action is 
a quite important topic, because one is trying to optimize the speed of the computation 
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and the width of the scaling region. In the lattice 1/N expansion, which we are discussing, 
another relevant criterion of choice is the possibility of performing analytic calculations as 
far as possible, in order to keep under control and possibly minimize the number of the 
altogether unavoidable numerical integrations. 

In a formal approach (and for one-coupling models), one might follow the original sugges- 
tion by Stone ||121|| , and adopt a lattice action defined by the kernel of the heat equation on 
the manifold where the fundamental degrees of freedom are defined. This choice has the ad- 
vantage of corresponding to the fixed-point (continuum limit) action in the exactly-solvable 
one-dimensional case [122]. However for our purposes it is definitely more convenient to 
consider actions that are polynomial in the fields. 

In passing we mention that a lattice Hamiltonian formulation was introduced by Hamer, 
Kogut, and Susskind ||123|| , and studied in the large- N limit in Refs. [124- 126 |. 1/N correc- 
tion have not however been evaluated. 

We shall discuss a number of possibilities, and finally focus on those formulations that 
meet the abovementioned criterion. 



A. Nearest-neighbor quadratic actions 



The lattice counterpart of Eq. ( |3.1| ) satisfying the request that only nearest-neighbor 
interactions be involved and no term higher than quadratic in any given field be present is 



where the iV-component complex field z n satisfies the constraint 



(6.2a) 



and we have introduced explicitly in the action a U(l) gauge field A ni(U satisfying 

Xn lf iK,n = 1 • (6.2b) 

This form of the action was introduced in the literature and its large N features were 
discussed in detail in Ref. [[/J]. Recently the case p v = (CP^ 1 models) of Eq. Q has 
been the starting point of many numerical simulations at small and intermediate values of 



N [127-130 



B. Nearest-neighbor quartic actions 



It is possible to write down a gauge-invariant lattice action without explicitly introducing 



a U(l) gauge field |5M52 



Sg4 — N[3 g ^ ^ [l — \z n+fJ- Z n \ ] 



(6.3) 
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This used to be a favorite version of lattice CP N ~ l models, especially because of the property 
that the N = 2 action is completely equivalent to the popular standard lattice action of the 
SU(2) w 0(3) nonlinear a model. is sometimes referred to as the "adjoint" form 

of lattice CP^ -1 models, and is characterized by possessing a large- N first-order phase 
transition, that does not however show up for any finite value of N. 

The mixed action obtained by combining (for j3 s = 0) and S^J is also solvable in 
the large- N limit |7[. More generally, the qualitative features of the phase diagrams for 
mixed actions can be explored, for finite N and in the 1/N expansion, by the mean- field 
technique | |133| - |135| ], which is an increasingly accurate description of the lattice models when 



one considers higher and higher space dimensionality. 

However large-iV studies and numerical simulations have shown that the approach to 
scaling of is very slow even for quite large correlation lengths, and the situation is even 
worse for what concerns asymptotic scaling, since the large- N (3 function gets contributions 
to all loops, in contrast with the continuum and the action ( p.l| ) (cf. Ref. ||134|| ). We shall 



therefore avoid any further effort concerning Eq. ( |6.3|) and its 1/N expansion. 

C. Next-to-nearest-neighbor actions 

It is obviously possible to formulate the models in terms of next-to-nearest-neighbor 
interactions: 

g(2) _ iV^^ j^/3 v [2 — Z n+ 2fj,Z n — Z n Z n+ 2f//\ 

n 



n,ii 



and similarly 



S$ = lNp e y %2[l-\z n+2(l z n \ 2 ]. (6.5) 



n,fx 



More generally, any combination in the form 

Cl S w + c 2 S {2 \ ci + c 2 = 1, ci > (6.6) 

is a lattice action belonging to the same universality class. Amongst these combinations, a 
special role is played by the choice 

£Sym = 4^(1) _ 1^(2) ? 

corresponding to the so-called Symanzik tree- improved version of the model ||136| , |137|| . The 
choice S Sjm corresponds to forcing a short- distance (ultraviolet) behavior markedly more 
similar to the behavior of the continuum action. As a byproduct, the leading logarithmic 
dependence on the correlation length of the deviations from scaling is removed and the 
scaling region is therefore enlarged. This makes Eq. (677) a natural candidate for numerical 
simulations of the models. 

In our discussion of the 1/N expansion, we shall try to carry the analysis of the general 
case CiS^ + c 2 S^ as far as possible, and we shall concentrate on when its peculiar 
analytic properties will become crucial to the development of our study. 
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7. 1/N EXPANSION ON THE LATTICE: EFFECTIVE ACTION, 
PROPAGATORS AND VERTICES 



In the context of the 1/N expansion, the crucial property that selects quadratic actions 
is the possibility of performing the exact Gaussian integration over the z fields at the only 
price of introducing a scalar Lagrange multiplier for the constraint 

Let us introduce the matrices 



4<^g ^ ^ (25m, n -^m,/i-^m+At,A*5m+2/^,n ^m,/i^m-/i,/i^m,n+2/i) 



(7.1a) 



(7.1b) 



and obtain the effective action in the form 



S cS = iVTrln 



ci-M m ^ + c 2 M$ n + — «n5m,r 



N 



(7.2) 



In order to perform the expansion, we must solve also the constraint involving the A„ )M fields, 
which can be done by setting 



A nijU = exp (itf ) 



(7.3) 



where 9' n „ is an (unconstrained) real field essentially playing the role of 9' in the continuum 
version. 

Assuming translation invariance of the classical fields a n and 9', we can write down 
the effective action in momentum space: 
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2/ 



A solution of the saddle-point equation is easily found to be 

9^ = 0, a = -im 2 , 
where is defined by the mass gap equation 
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d 2 p 



1 



2f J (2tt) 2 ci £ „ 4 sin 2 ±p M + c 2 V sin 2 Vil + m 2 7 (2vr) 2 p 2 + m 2 



d 2 p 
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(7.4) 



(7.5) 



(7.6) 



31 



we adopted the standard notations 

p M = 2sin|?v> ft = ^P% P 4 = ^pJ, (7.7a) 

P 2 = ci ^ 4 sin2 f?V + c 2 X] sin2 Pv = J ^2pI ( c i + c 2 cos2 §2V) = P 2 ~ i c 2P 4 , (7.7b) 

(7.7c) 



2„ 



d 2 P = r dpi r dpa 

(2tt) 2 - J_„ 27T 7_ w 27T ' 

Going back to coordinate space, taking the second functional derivatives with respect 
to the fields a n and 9' , and evaluating them at the saddle point defined by Eq. (|7.5|), we 
obtain a representation of the lattice propagators of the effective fields: 

a mW = fw^ =J =i ' ( 7 - 8a ) 

J p +±k + m\p-\k +ml 

! /" d 2 p 2C1COSP,, + c 2 cos 2p M (l + cos k^) 

<W*) = (i + x /} ^ y _ WT ^ I 

(2ci sinp^ + c 2 sin 2p^ cos \ k^ (2ci sinpj, + c 2 sin 2^ cos 



p+\k + m 2 ^ fp— + ml 



(7.8b) 

+ m 2 J 

where the 9' field has been rescaled to 9^ = 9' /(l + xf). 

The effective vertices are obtained by taking higher order functional derivatives of the 
lattice effective action. They are in the form of one-loop integrals of the fundamental fields 
propagators. Computations can be cumbersome, but they are conceptually straightforward. 
The Feynman rules needed for the computations are summarized in Fig. [12|. We used the 
notation 

^3|t(p, k) = 2ci sinp^ + c 2 sin 2p^ cos (7.9a) 
Vfyip, k) = ci cos^ + c 2 cos 2p^ cos 2 \k^. (7.9b) 

As discussed in Sect. ^|, the vertices are gauge-invariant for all values of x. As a consequence, 
the following relationships hold: 

2 2 

2JVMp,*0 =P+\k -p-\k , (7.10a) 
m 

2MV(P, fc) = ^(p+|fc, fc) - ^(p-ifc, fc). (7.10b) 

The notations we have introduced allow the following representation of the effective 
vector propagator 

d 2 p V ifl (p,k) 
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(27r) 2 p 2 + m 



-■.^-gi/g yy?' ( ,u) 
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By applying Eqs. ( |7.10| ), we may check that 

(27r) 2 p 2 + TOq 



^exp {-\ik^) hpi^lS^ J 



exp [±i(kp - K)\ I — — — * 7=* — ) = 0, (7.12) 

l Z7r J (p 2 + m 2 ,) (p+A: +wig 



as an effect of gauge invariance: when k = the inverse propagator zLgJ (A;) is transverse 
on the lattice, i.e. it vanishes when contracted with e~ lk ^ — 1 or e %kv — 1. It is therefore 
convenient to restate Eq. (|7.11 ) in a form reminiscent of Eq. (|3.9| b): 



A TLM = *mM + A mW C(k), (7.13) 



where 



W = <W - exp - M] (7.14) 

k z 



(2?r)2 - 2 + m 2 > i 7 - 15a ) 



enjoying the important property 

A-](0) = 0. (7.16) 

It is worth noticing that 

independently of and fi. However even in this simplified case one cannot trivially identify 
x with its continuum counterpart, even in the scaling regime, and, for finite values of /, 
xi att must be tuned in order to recover an assigned value of x cont . 



8. INTEGRAL REPRESENTATIONS OF LATTICE PROPAGATORS 

We mentioned in Sect. ^ that some powerful analytical techniques can be applied in some 
specific version of the model. In this Section we introduce the first of these techniques, the 
use of integral representations of lattice propagators. We could only apply this technique to 
the case c 2 = (nearest-neighbor interactions) which we shall discuss in great detail. 
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When C2 = dramatic simplifications occur already in Eqs. ( fT.8|) . By simple manipula- 
tions and use of the gap equation 
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(8.1) 

(8.2a) 
(8.2b) 



The first analytical result concerns Eq. (|8.1| ), that can be cast into the form 
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2tt 1 + mg/4 \l + ml/A 



(8.3) 



where K is the complete elliptic integral of the first kind. It is also possible to evaluate in 
closed form the inverse propagators along the principal diagonal of the momentum lattice, 
i.e. when k\ = k 2 = I. We change variables to q\ = \{j>\ + P2) and q 2 = \ij>i — P2), and 
notice that, for any periodic function f{p\,P2) we have 

In the case fci = fc 2 = I in particular 

j3 2 = 4(1 — cos qi cos g 2 ), 
' — 2 

p+k = 4(1 — cos(gi + I) cos (72). 



(8.4) 

(8.5) 
(8.6) 

Therefore the g-integrations of Eqs. (|8.2| ) can be performed in terms of the standard complete 
elliptic integrals K, E, and II (we use the conventions of Ref. ||138|1 ); the result is 
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(8.7a) 



1 E 



1 + m 2 /4 



+ 1- 



COS 2 \l J 7T 



-K 



1 + m 2 /4 



+ 



cos 2 hi 



n 



cos 2 \l 



(l + m 2 /4)V n V(l + m 2 /4) 2 ' l + m 2 /4 



(8.7b) 



In the general case, we had to resort to an integral representation of the inverse propa- 
gators. Let us make use of the standard Feynman parameters to write 



1 



1 



p 2 + m\ p+k + itLq 



dx ■ 



1 



p 2 {l — x) + p+k x + m\ 



1 2 



(8.8) 
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and notice that, via trigonometric identities, 

(Pn+K) 2 =Pl{ 1 - \K) +K + 2 sin ?v sin K ■ 

Let us now change variables to g M = p^ + z^, where z^ is defined by the relationships 

x sin k u 



(8.9) 



sin z,, 



1 2 H 



1 — x(l — x)k 2 
Eq. ( |3.10D implies the identities 



1 — x(l — x)k 



pl(l -x) + {p,+k,fx = fiy/l-x(l- x)kl + 2 - 2y/l - x(l - x)kl 
sin g M cos ^k^ + (1 — 2x) cos g M sin ^k^ 



sin (p M + \k^) 



1 — x(l — x)k 2 



Substituting Eqs. (p>.8|) and (|8.11|) in the relevant one-loop integrals we obtain 



dx 



d 2 q 



1 



A 



^ (27T 

! _ 4 + ml 



(0) 



- 1 - / da; 
2/ Jo 



d 2 g c - J] 6 M sin 2 g M 



2 • 



(2tt) 



1 2 ' 



where in order to simplify the notation we have defined the auxiliary variables 



a M = 2y 1 - x(l - x)k*, 

(1 - 2a;)4 2 + kl - 4 

0/i = x 

1 — x(l — x)A; 2 



We can now exploit the relationships 



2 d 2(1 -2s) 
b = a 



5.10) 



5.11a) 
5.11b) 



5.12a) 
5.12b) 



(8.13a) 
(8.13b) 

(8.13c) 



(8.14a) 
(8.14b) 



to perform an integration by parts in the variable x in Eq. ( |8.12| b) and obtain the more 
convenient representation 
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d 2 q 



4 + ml - Y^fj, a M cos 



_LA I d2q 



aj sin 2 



(8.15) 



Trivial algebraic manipulations (involving repeated use of integration by parts in the mo- 
mentum variables) lead finally to the form 



- 2 f AX (l^l 

Jo a\ci2 



d 2 q 



4 cos 9i cos q2 



1 2 " 



(8.16) 



Starting from Eq. ( |8.12| a) and ( |8.16| ) we can now perform the momentum integrations. One 
momentum component can be integrated easily, thanks to the relationship 



* de 



i 



_„. 27T 6 - acosfl y/b 2 - a 2 
and its parametric derivatives. As a consequence we obtain 

" n dqi 4 + ml — a\ cos qi 



(8.17) 



(a) 



dx 



o 



2vr 



(4 + ml — ai cos gi) 2 — ai 



3/2 



Jo ai J-n 2 ^ 



\0) 



4 cos qi 



(4 + ml — ai cos qi) 



3/2 



(8.18a) 
(8.18b) 



The q\ integrations in Eqs. (|8 . 1 8 ) are easily reducible to standard elliptic integrals [138 



Without belaboring on the straightforward algebraic tricks involved in the derivation, we 
present the final result in the following form: 



, _ 4 + A 1 1 C 3 £(C) 

{a) " 4tt J aX {a x a 2 y/ 2 1 - C 2 ' 



where 



p 4 + m| 



C 



dx 



;i - 2xf 

{a ia2 fl 2 



2 IV 3 



C 3 £(C) 



+ 2C(E(C)-K(Q) 



(8.19a) 
(8.19b) 



ACL1CL2 



(4 + m 2 ) 2 -(a 1 -a 2 ) 2 ' 



(8.20) 



Eqs. Q8.19Q lead to an enormous computational gain: from the numerical point of view, the 
values of the elliptic integrals can be routinely generated with high accuracy, and therefore 
the number of numerical integrations is reduced from two to one; moreover, the x integration 



36 



is much more regular then the original momentum integration: if we exploit the explicit ex- 
pressions of A^ and A^ on the principal diagonal Q8.71 ) to further regularize the integrand, 
a moderate-size Gauss-Legendre integration suffices to produce very accurate results. From 
the point of view of analytic manipulations, we can exploit the knowledge of the asymptotic 
expansions of the elliptic integrals to simplify dramatically the otherwise quite complicated 
problem of generating an asymptotic expansion of the propagators in the small mo (scaling) 
region. This expansion in turn is the essential ingredient in quantitative computations for 
lattice models in the scaling (i.e. field-theoretical) regime. We shall discuss this point in the 
next Section. Let us only recall here the relevant expansion formulae of the complete elliptic 
integrals, in the regime k ~ 1, k! — \/\ — k 2 <C 1: 



We are interested in evaluating physical quantities in a lattice model possessing a non- 
trivial (continuum) field theory limit when the ultraviolet regulator (1/a, where a is the 
lattice spacing) is sent to infinity while properly tuning the lattice coupling / according to 
the renormalization-group equation of the model. In such a model, lattice expectation val- 
ues for finite a and f(a) will necessarily receive contributions from the irrelevant operators 
included in the lattice definition of a physical quantity. In order to isolate the contribution 
that will survive in the continuum limit, i.e. the scaling part of the expectation value, we 
must be able to perform an expansion in the powers of the lattice spacing. Within the 1/N 
expansion of lattice spin models, there is another dimensionful parameter that can be em- 
ployed to classify the relevance of the contributions to any given expectation value: the bare 
large- iV vacuum expectation value ( "large- N mass"). Since the continuum limit is the 
limit of infinite correlation length, it is also the limit of vanishing m^: contributions carrying 
higher powers of with respect to the scaling contribution will be irrelevant. Actually, 
since the physical value of the mass in these asymptotically-free models is not strictly zero, 
we must be careful not to set = and to get our results in the form of asymptotic 
expansions in which 1 ^> |l/ln(mQa 2 )| ^> m\a 2 . 

The detailed analysis of the procedures for the evaluation of physical quantities in the 
scaling region, and the related problem of regularization of the infrared divergences generated 
by the naive — > limit, will be discussed in Sect. O. Preliminary to such a discussion 
is however the asymptotic expansion for small values of of the basic ingredients in the 
evaluation of physical quantities, i.e. the effective propagators themselves. Due to their 





(8.21a) 




9. ASYMPTOTIC EXPANSIONS IN THE SCALING REGION 
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original definition as one-loop integrals over Feynman propagators, the inverse effective 
propagators Z\7% and ^7g^ u will have formal asymptotic expansions in the form 

oo 

[MQ + B n {k) lnmjj] {mlY . (9.1) 

n=0 

We may however recognize from Eq. ( |7.6|) that the coupling constant itself admits such an 
asymptotic expansion in the form 

oo 

P = 27 = E K + & lnm o) KT , (9-2) 
and Eq. ( |9.2| ) can in principle be inverted to 

oo 

\nml = Y,(ln + S n (3){m 2 ) n . (9.3) 



n=0 



As a consequence, a perfectly acceptable form of the asymptotic expansion of the inverse 
propagators is: 

oo 

= E W a) (*) + M Q) (^)] W) n , (9.4a) 

n=0 

oo 

# = E + /^f M)" ■ (9.4b) 

n=0 

The expansion ( |9.4| ), in comparison with ( |9.1| ), involves no extra effort, but it shows the 
great advantage of expressing final results in a form that makes direct contact with the 
lattice weak-coupling expansion of the models, and allows explicit checks of commutativity 
between weak-coupling and 1/N expansion. 

We obviously aim at a systematic way of generating the coefficients A n and B n . We shall 
discuss here a general technique that does not rely upon any specific feature of the one-loop 
integrals involved, and in the next Section we shall present a different approach specific to 
the case discussed in Sect, p] when an integral representation is available. 



A. General technique 

The general technique is worth a detailed analysis, because it will be very useful also in 
the evaluation of expectation values of physical quantities. It will also allow us to show the 
connection between sharp-momentum regularization and lattice formulation. It is therefore 
convenient for our purposes to discuss the case of a one-loop lattice integral with arbitrary 
propagators and vertices. The general form of the lattice integral is 

r d d p 

I(k;m a) = J j—^ F(k;m a,p), (9.5) 
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where k is the external momentum (or a collection of external momenta) and we use the 
notation 

6 d d p d f b dpi f b d d p d f b dpi d f a dpi 

(2^=ny_ 6 ^' i j^=uj_ b ^-iij_ a ^- ( g - 6 ) 

When comparing with the continuum formulation, we must keep in mind that powers 
of a have been included in the definitions of F and / in order to make them dimensionless. 
Formally we can rewrite Eq. ( |9.5| ) in the form 

f n d d p f 00 d d p 

I(k;m a)= ——- F(k;m a,p) + a d — — ; F(k;m a,pa) 
Jo i 2 ^) Jo ( 27r ) 

r°° d d r) 

~ a j^TFfam^pa), (9.7) 
Jo { Z7l > 

where TF is the Taylor series expansion of the function F in powers of mod around m^a = 0. 
In order to turn Eq. ( |9.7|) into a mathematically rigorous statement, we must show that, 
order by order in a, all ultraviolet and infrared divergences are explicitly canceled. Let us 
therefore introduce the truncated Taylor series expansion T j (including all the powers of mo 
up to 771q) and notice that 

f n d d p f°° d d p f°° d d p 

/ 7T~u F ( k ; m oa,p) + a d / 777-^ F(k; m a,pa) - a d / — -- TjF(k; m a, pa) 
Jo \ Z7T ) Jo [ Z7T ) Jo \ Z7T ) 

r d d p f n ^ a d d p 

Jo (2tt) Jo (2vr) d 

+ a d ——^F(t,moa,pa) + a d / -— - (1 - Tj) F(k; m a,pa) 

Jo Jir/a 

r d d p f°° d d p 

= jTr-Td ( X - Tj)F(fc;m a,p) + a d / — — ^ (1 - Tj) F(k; m a,pa) + I(k; m a) 

Jo K^) Jir/a \ Zl1 ) 

= I(k;m a) + 0(a J+1 ). (9.8) 

Indeed the first contribution to Eq. ( |9.8| ) is trivially 0(a J+1 ), and the same can be shown 
to hold for the second contribution by noticing that 

F(k;m a,pa) = aPFj(k; mp,p), 

j 

Fj(k- Xmo, Xp) = X^F 3 (k- m ,p); (9.9) 

therefore we have 

a d ^- d (l- Tj ) Fj (k;mo,p) = 0(a J+1 ). (9.10) 

J-w/a 

Eq. ( |9.8|) paves the way to two major strategies in the regularization of lattice integrals. 
The first strategy, till now most popular, stems from the observation that the last term in 
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the l.h.s. of Eq. ( |9.8| ) vanishes in dimensional regularization. One may therefore adopt the 
definition 



I(k: nina) = lim 

d^2 



d d p 
(2ir) a 



F(k; m a,p) + a° 



d d p 
j2nY 



F(k; m a,pa) 



(9.11) 



When expanding in powers of a, the first contribution corresponds to infrared-divergent 
massless lattice integrals and the second to ultraviolet-divergent massive continuum inte- 
grals; in dimensional regularization infrared and ultraviolet poles cancel exactly and we are 
left with a finite result. However one may also notice that the last term in the l.h.s. of 
Eq. (|9.8| ) is reminiscent of the structure of counterterms in the sharp-momentum regulariza- 
tion scheme. More precisely, we can decompose the expansion TjF into contributions which 
are ultraviolet divergent (relevant), contributions which are infrared divergent (irrelevant), 
and contributions which are ultraviolet and infrared divergent (marginal). Symbolically we 
may write 



rp(UV) + T (iR) + T (o) 



(9.12) 



We can introduce an arbitrary cutoff M 2 and split the counterterm Tj ^ into two separate 
integrals. Finally we can express the original integral in terms of two separate contributions 
whose singularities are now independently regularized: 



I(k; m a) 



o (271T Jo 47r 







-j— T^?F(k; m a, pa) 

4:71 



+ 



F(k; mod, pa) 



— a 



Tj UV ' ) F(fc; m a,pa) — a 2 

4li 



o 

oo 
M 2 



50 ay 

4tt 

^-T5 )) F(fc;m a,pa) 

47T 



(9.13) 



Eq. ( |9.13|) is the starting point for the series expansion in powers of a. 

The first contribution in brackets is ultraviolet-regular. The lattice integral and the con- 
tinuum integrals are separately infrared-singular, the infrared divergences canceling out in 
the sum. In analogy with the standard notation for ultraviolet regularization (of. Eq. (|5.1|) ), 
we always assume that such combinations of (massless) separately infrared-divergent inte- 
grals stand for e.g. 



d 2 p 



'latt 



(P) 



d 2 p f dp 2 

Olatt(p) - 0cont(p)) - / 0cont(p), 



(27T 



■R 



4ix 



(9.14) 



where 1Z is the region comprised between the circle of radius A and the square of side 2tt 
(assuming A > rr^/2). The integral of the difference is infrared-regular by construction. 

The second contribution in brackets is infrared-regular and ultraviolet singularities are 
removed according to the prescriptions of the sharp-momentum scheme, which means that 
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the marginal (scaling) component of I(k; mod) receives a contribution exactly equal to the 
one obtained in the corresponding continuum model: all lattice effects are included in the 
first term. Moreover, as we shall immediately show, by applying Eq. ( |9.13| ) directly to the 
gap equation it is possible to fix the value of M 2 in such a way that the sharp-momentum 
coupling constant can be identified (for large N) with the lattice coupling. 



B. Expansion of the gap equation 



For our purposes it will be sufficient to truncate the expansion in powers of a to the 
second nontrivial order, i.e. J = 2. The regularized form of Eq. (|7.6|) is 



d 2 p 



2f J (2tt) 2 m 2 a 2 + V 4 sin 2 ±p M - c 2 £4 sin 4 ±p M 



+ a z 



dv 2 1 

~g h _ <-£\ 

47r m\a 2 + 4 sin2 \ a Vv - C2 4 sin4 \ a Pv 



and, expanding to 0(a 2 ), we obtain 



1 



+ 



* d 2 p 
dp 2 

4:71 

dp 2 

47T 



2™2 



am 



o 



1 a 2 (c 2 +|)£ p 



+ 



4 



2^2 



p 2 + m 2 4 (pa + ml) 

1 m 2 + 1 a 2 (c 2 + i) p£ 1 a 2 ™ 2 , (c 2 + |) £„ P', 



p^ 



(p 2 



4 



(p 2 Y 



,2\ 3 



According to Eq. fl9.13|) we can write 



* d 2 p 
2f ' J (2tt) 2 



a 2 ml 



P 2 ~ ic 2 p 4 



+ 



dp 2 

47T 

00 dp 2 

47T 

°° dp 2 r 
47T 



(p 2 - ±c 2 p 4 ) 
1 1 a 2 m 2 . ^ E M P 

l C 2 + 3 J 



+ 



dp 2 wiq 
4^ (p^f 



2 p 2 



(p 2 



3 a 2 (c 2 + |) (p 2 ) 2 3 



+ T77 



p 2 + m 2 ' 16 (p 2 + m 2 ) 2 16 ° ^ + ^ 



M 2 



3 a 2 mQ 



(<» + t) 



p^ 8 p^ 

The continuum integrals can be explicitly evaluated and the final result is 
" ~ d 2 p 1 /"^dp 2 1_" 

47T p 2 



1 

2f 



o (2vr) 2 p 2 - I C2 p 4 



1 , M 2 

47T 







(9.15) 



(9.16) 



(9.17) 
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— am 



d 2 p 



+ fl2m o ^ ( c 2 + §) 



1 - 2 In 



M 5 



dp 2 1 

47T (p 2 



dp 2 1 , E„pj 



4vr 2 



(9.18) 



Eq. Q9.18 ) deserves some comments: 

A rescaling p — > p/a in the terms originated by Tj must be performed after the 
expansion in powers of a has been accomplished. It may sound arbitrary, because we are 
working with divergent quantities without an intrinsic scale, but it is mathematically sound 
as one may recognize by going back to the proof of regularization. 

The effect of a Symanzik improvement is made apparent by the cancellation of the 
logarithmic dependence on m in the first irrelevant contribution when we choose the value 
c 2 = -§. 

In order to identify sharp-momentum and lattice couplings we must choose M\ such that 



d 2 p 
(2tt) 2 



Ml 



Att p 2 



(9.19) 



implying 



1 1 M? 
— = —ln—k 



2f 4vr 



m 



i 2 2 

+ a m Q 



o 



ai{c 2 



l + 3c 2l M\ 



32tt 



In 



+ 0(m 4 c 



(9.20) 



The choice implied by Eq. ( |9.19| ) will allow us to express our results in a form immedi- 
ately comparable to the lattice weak coupling expansion, because the difference between 
the couplings is purely nonperturbative, as shown by Eq. ( |9.2(J| ). Any other choice would 
correspond to a finite renormalization and would require a perturbative adjustment of / in 
order to recover standard perturbation theory. M\ and a.\ are plotted as functions of c 2 in 



Fig. 13. Let us notice the special values 



^lI C2= o 



32, 



M 2 1 , /„ 

h\c%=— 1/3 



c 2 =0 



32tt 



c 2 =-l/3 



= 17.68967299, 
-0.00479767. 



(9.21a) 
(9.21b) 



C. Expansion of the propagators 

In practice, when willing to compute the expansion indicated in Eqs. ( |9.4| ), much work 
can be saved by the following observations: 

The functions B n (k) are related in a simple way to the coefficients of the 1/e poles 
resulting from the continuum integration in the representation fl9.ll ). These coefficients are 



related to essentially trivial one-loop continuum integrals that can all be computed in closed 
form. Therefore the functions B n (k) can be expressed in terms of elementary functions of k 
and do not require an integral representation. 
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Moreover, let us assume that the integrand in the representation ( |9.5| ) of an inverse 
propagator A' 1 be a function F^(k; m ,p). By shifting the integration variable it is possible 
to symmetrize F to F s J m {k\ mo,p), where F sjm is an even function of with all singularities 
located in p^ = ±^& M . Let us now consider the function 



d 2 p 
(2^)2 

d 2 p 
(2^)2 



Fr(k;m ,p)-pY, B n(k) (mjj) B 



n=0 



F s J m (k;m ,p) - - 



+ 



p+\k p-\k 



n=0 



(9.22) 



By construction, according to Eq. Q9.4p , this function is analytic in m around m 
Therefore we immediately obtain the following integral representation of A n (k): 



Mk) = 1 / 



d 2 p d n 



(2vr 



I 2 d(mg) n 



F s J m {k;m ,p) 



+ 



p+^/c + m 2 , p— | A; + 



z=o 



(9.23) 



mg=0 



The symmetrization procedure insures us that all the singularities in the integrand are locally 
canceled out. We are still left with singularities regularized by counterterms of the form 

cos 416 

(9.24) 



-,2r 



where q and 9 are the polar coordinates in the p plane with center and I 7^ 0; these 

counterterms are implicitly introduced whenever needed, or equivalently polar integration in 
two small circles around p^ = is understood performed first. In practice, we shall only 

be interested in the functions A (k) and Ai(k), whose integral representations we explicitly 
quote: 



A (k) 



A^k) 



d 2 p 

d 2 p 
(2^ 



F s r(k;0, P )-- 



dF s J m (k;m ,p) 



p+\k 



p-\k 



B {k) 



(9.25a) 



dm 2 . 



+ 



mj=0 



p+\k 



■2\2 



+ 



p-\k 



2\2 



B (k) 



+ 



1 



p+lk p—\k 



Bi(k) 



(9.25b) 



Hence for our present purposes we shall only have to work out the functions Bo(k) and 
Bi(k). Let us illustrate the procedure by evaluating B^\k) and B[ a \k). Our starting 
point will be the expansion in powers of a of the continuum integral 
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x 



fl2 + 4 Sin2 _ C2 X] 4 Sin4 (1°^) 

A 4 A 4 

o fl2 + 4 sin2 i^ ap v + k ^) ~ ° 2 ^2 4 sin4 (s ^ + M 



(9.26) 



where the factor 2 has been included in order to take into account the effects of the (excluded) 
expansion around the second singularity in p = —k. Since we are interested in the expansion 
to second nontrivial order we may replace Eq. (|9.26|) with 



2a £ 



oo A d 



dp 



1 fl2 S M v\ (ci cos kf, + c 2 cos 2kp) + ml 



o (27r) d pa 2 + m^a 2 I /c 2 



+ 



( 2c i sin fc M + C 2 sin 2fc P 



(9.27) 



Within the desired approximation however 

roc A d„ n „ X /-oc d d p 



o (2vr) d p 2 + mg d J (2n) d p 2 + m 2 , 



_m 2 /"°° d d p 



o (27r) d p 2 + m 2 , 



(9.28) 



Moreover, we are only interested in the pole part, and we can therefore replace Eq. (p. 27 ) 
with 



oo A d 



dp 



o {2n) d p 2 + m 2 I k 2 



1 

— — m^a 



2„2 



1 _ \ J2n ( c i cos + °2 cos 2fc p 



(k 2 y 



i Z) M ( 2c i sin fc M + c 2 sin 2 K 
2 



(9.29) 



However, 



oo id 



dp 



o {27r) d p 2 + m 2 



is an exact representation of the pole part in the asymptotic expansion of (3, and therefore 
we can read B ( a) and B[ a) directly off Eq. ( ggg ) and find 







2 1, 



d(«) sin 2 (ci + 4c 2 cos 2 i/c^) ^ sin 2 fc M (ci + c 2 cos 2 
^i = ~ 2 Z^ 77772 4 2^ 



An important representation of Eq. (|9.30|b) is 



[k 2 y 



(9.30a) 



(9.30b) 
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-°1 — -°o 



1 Q Q T 72 , Sin2 \ k V ( C l + 4C 2 COs2 \K) 2 2 

j£Wn* 2 + £ fc 2 — ~F 



(9.31) 



The computation of _Bq and is in no sense conceptually more involved than the 
computation we have just presented. However many more terms have to be taken into 
account and a few algebraic tricks exploiting explicitly the fact that there are only two 
vector components in two dimensions have to be employed in order to simplify the result. 
Let us only quote the final result: 



J BS e) = 2^(c 1 + c 2 cos 2 |^)-2^ 



kl ( c i + c 2 cos 2 -hp)' 



2 n ( Ci + ° 2 c ° s2 



k 2 



C\ + 4c 2 cos 8 cos (ci + c 2 cos k 



k 2 



To\2 



(k 2 ) 



J2 U sin 2 \k v (ci + c 2 cos 2 \k v ) 



c\ + c 2 cos 2 \k 



- sin 2 \k^ 



IE 



sin 



2 1 



kf, (cj - 4c 2 cos 2 |A^) 2 



A' 2 



Ci + c 2 cos 2 1^ 



21 8ci sin 2 
ci + 4c 2 cos 2 ifc 



A' 2 



(9.32a) 



(9.32b) 



Finally we mention that it is possible to define an asymptotic expansion also for the 
functions A7 e ~, {k) defined in Sect. |7|. The form of the expansion is 



n=0 



and it is easy to obtain 
D 0fJ , = ci + c 2 cos 2 



a 



d 2 p cip 2 + 4c 2 ^2 v pl cos 2 t}p„ cos 2 \k^ 



(2tt) 2 

-Dim = I ( c i + 4c 2 c °s 2 \k^) , 



C 



d 2 p 



ci + 4c 2 cos 2 \k» cip 2 + 4c 2 J2u Pi cos 2 \p v cos 2 



.2 1. 



p 2 



{p 2 



(9.33) 

(9.34a) 
(9.34b) 
(9.34c) 



(9.34d) 

All these results take a particularly simple form in the case c 2 = 0; in this case we obtain 
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2 




-(- 

V 


Bj°\k) 
Bf\k) 


= 2, 

1 

~ ~2 




= 1, 




4 k 4 



(9.35a) 

(k) (~£W^ 2 + \~ , (9.35b) 

(9.35c) 

fiE 5 Alnfc 2 + U^ ) , (9.35d) 



p (p)^ \ 2 4 ' ^ 



C 0/1 = -J, ^im = I 5 C 1;J = 0, (9.35e) 



and one can easily construct the representations of Ao(k) and Ai(fc) by a trivial application 
of Eqs. ( |9.25| ). We must observe that these representations are well defined and they can 
be used to compute numerically A (k) and Ai(k); however, the integrands are plagued with 
a very singular behavior, especially for small k, and it is very hard to perform accurate 
numerical integrations. 



10. ASYMPTOTIC EXPANSIONS USING INTEGRAL REPRESENTATIONS 

The methods developed in the previous Section apply perfectly well to the case C2 = 
0. However this is not the most convenient way of performing the asymptotic expansion 
when we possess an integral representation like Eqs. (|8.19|) . In this case we start from the 
observation that 

I(k;m a)= / dx G{k^ m a, x). (10.1) 
Jo 

Moreover, if we perform an homogeneous expansion of G in the powers of k^ and mo, the 
x-integration of the resulting coefficients can be performed explicitly. As a consequence, we 
may write 

I(k;moa)= / da; 1 — Tj G^h^a, m a, x) + / dx T^G^^a, m a, x), (10.2) 
Jo L J J 

where the unconventional notation k„a indicates that the parameters of the Taylor expansion 



are k^a, not k^a. The first integral is now regular up to 0(a J ), and the integrand can 
therefore be expanded in powers of tuq. The second integral in turn can be performed 
analytically and the result, which is a nonanalytic function of mo, can be expanded in an 
asymptotic series in the powers of mo- This procedure is rather cumbersome, but the payoffs 
are very high, as we shall see in the following. 

Let us present the essential ingredients of the computation and sketch the derivation of 
the asymptotic expansion in the case of Let us define 

(£ = £js(i-s), g 2 = E* ? 4 = E<' ( 10 - 3 ) 
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and expand all relevant quantities in powers of wig and A; 2 : 



c 2 



1 - \ (m 2 + q 2 ) , 
1 _ m 2 + q 2 / m 2 + q 2 ^ 



K(C) « — In 
vs; 2 32 



32 

1 . ml + q 2 ml + q 2 ^ ml + q 2 ^ml~ q 2 



16 



32 



16 



8 q 2 + ml 



(10.4a) 
(10.4b) 

(10.4c) 
(10.4d) 



The relevant integrals are: 



dx 2 , f + 1 

— 7T = -TT^: In , 

q 2 + m 2 k 2 ^ £-1 
dxln^±^ = |lnC±i-2-ln^, 

32 e-i m o 



dx ■ 



o (g 2 + mg) 2 (A; 2 ) 



1 



1 3 

— + 7T 

2^2 2 



(10.5a) 
(10.5b) 
(10.5c) 



where 



1 + 



Ami 
k 2 



(10.6) 



Then, starting from Eq. ( |3.19| a), we obtain 



1 

^ JO 

1 

2^ 



dx 



m 2 , + q 2 



ml + q 2 q 2 + ml 



32 



4 4(ml + q 2 )_ 



1 ln^ + 1 



V^ 1 



1_2£ 

16 | ' 8 ' 16 m 2 



1 1, 32 
- + — In — * 



1 k' + i+1 3 | 2 - 1 

■ In 1 



(10.7) 



8 (P) 2 V K* < I ^ 2 

Therefore by expanding both sides of Eq. (|10.7| ) we can introduce the regulator in the form 



+ 



1 

2^ 

m 



, 1 1 i fc „ 1 m l 
dx — + — — In — + 2/5 — + — 

2g 2 2wk 2 32 /c 2 2tt 



dx 



1 3 1 



o 



2nk 2 



k 2 8 2 (p) 



A- 2 1 
, .In — + 2i 
2 ' 32 k 2 



1 1 A; 4 



1 A; 4 1 



16 q 2 4 (A^ 2 fl 2 



4 4 (fc 2 ) 



+ 



A' 2 



4 1 



(fc 2 )^ 



(10.8) 
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As a consequence, the known values of Bq and B± are reproduced, and we find the 
representations 



.4, 



(o.) 



A 



(o) 



+ 



+ 



— / (li- 
ft 



— I dx 

47T 
f 



Co 3 E(Co) 



1 , k 2 
in—, 



( ai a 2 ) 3 / 2 1 - Co 2 VJ 2ftfc 32 



Co 3 / e(Co) + 4c 2 r^(Co) - mco) 



(10.9a) 



Aft 
1 

27 



dx 
2 



(aia 2 ) 3 / 2 \ 1 — Co a i a 2 

1 3 1 1 q 4 

+ 



(g2)2 8 g 2 2 (g2)3 

1 1 k 4 1 ln k2 
J^f~^~4jk 2 Y F n 32 



2 3 
I 2 ~ 8 



2£(Co) 
(1 " Co 2 ) 2 



1 k 4 



2 (A; 2 )' 



(10.9b) 



where 



Co 



4aia 2 



16 - (ai - a 2 ) 



2 ' 



(10.10) 



By repeating the analysis for A,l we find that 



ZT 1 « — 
W 4vr 

1 

~ 27 



dx 



ein 



(1 - 2x) 
ml + g 2 

1+1 



2 r 



g 2 + ml 1 



^ 2 + — A: 2 (3-2f 2 ) fin 

e-i i6 1 



l-i 



32 
5 



4 ml + q 2 



JL^_1^2 ln 32 



24 



16 



1 fc 4 /(l-^)(l-5a i e+i , 13 



+ 8 jfc2 



4f 



In 



(10.11) 



and expanding both sides of Eq. fllO.lip we can introduce the regulator in the form 



irk 1 Ik 2 n mi 
— dx — + — In— + 2/3 + — ^ 
2ft J 2q 2 2ft 32 2tt 



dx 



A- 2 



+ 



3 k 2 Ik 4 



+ 



m5 
2tt 



1 k 4 



2 (g 2 ) 2 q 2 16 g 2 4 Pg 2 
I I 1 k 4 



+ m o 



4- 



A: 2 
I 1 



2 (P) 2 

£.4 



A: 2 

In — + 2 , 

32 A: 2 4 4 (A;2) 2 



A; 2 2 (jfc2y 

The values of Bq and B± are reproduced and we obtain the representations 



(10.12) 



-4, 



(0) 



4ft 2 f 1 ^ (1 - 2x) 2 
7T 7 (aia 2 ) 5/2 



Co 3 £(Co) 
1-Co 2 



+ 2Co(£(Co) - K(Co)) 
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- . k 2 1 , k 2 

— ax — - H in — 

* ' 2q 2 2vr 32 



(0) 



1 

TV 

4k 2 

71 



(10.13a) 



o 

1 (1 - 2xf 
dx- —Co 



dx 



( ai a 2 )V2 
-2xf 



3 +E{( )-2K(Co) 



'aia 2 ) 7/ 



2 



+ — I dx 

2tt 



1 

2tt 



r, 1 ~ Co + Co 771 / > 

1 3 1 



2-Co 2 
1-Co 2 



^(Co 



2 g2 p 16 q 



1 

2 " 4 




k 2 n 1 1 1 fc 
In — + 2 — + - - - ^r-^ 
32 k 2 4 4 (P) 2 



;i0.13b) 



The contour plots of y4 °^ and A Q a ' in the fc-plane are drawn in Fig. [14. 

We finally present the results of the evaluation of the inverse propagators along the 
principal diagonal of the momentum lattice, i.e. when k\ = k 2 = I. In this case direct 
manipulations of Eqs. (8/7) lead to 



4 a) (M) 
Af\l,l) 
Af{l,l) 
Af\l,l) 



1 



167r sin 2 






1 




647r sin 4 


¥ 


1 


1 




27 


cos 


¥ 




1 





87r sin 2 |Z 



l-cosi/ 

cos f:t m f- , 

2 1 + cos \l 

cos 2 i/ - (2 - cos 2 i/) cos U In ^ — C ° S 2 } 
2 2 2 1 + cos \l 

1 — cos \l 
In h-2 



q -| -i 1 COS r) W 

2 — COS 7:1 + COS tj/ In f- 

2 2 1 + COS iZ 



(10.14a) 
(10.14b) 
(10.14c) 
(10.14d) 



11. EVALUATION OF PHYSICAL QUANTITIES IN THE SCALING REGION 



The regularization ( |9.13|) of lattice integrals in the scaling region and the asymptotic 



expansion ( |9.4| ) of the effective propagators are the essential ingredients for the evaluation 
of the scaling contributions to any physical quantity to 0(1/N), involving one- loop integrals 
over the effective fields propagators. 

Actually there is no real need of evaluating separately every physical quantity: whenever 
two operators have the same (possibly anomalous) dimension, their ratio is a pure number 
that is scheme-independent, and therefore it can be computed in the simplest available 
scheme; only one of the two operators, possibly the easier to compute, must be evaluated 
on the lattice in the scaling region. These statements obviously rest on the existence of a 
renormalization-group behavior in the scaling regime of the model involved, but they can 
be explicitly verified in the 1/N expansion, as we shall see in the following. 
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A. Lattice correlation length 



A prototype dimensionful quantity is the inverse correlation length (mass gap), that sets 
the scale for all quantities having canonical dimension. For sufficiently large values of x, as 
discussed in Sect. [| it is meaningful to extract the mass gap from the asymptotic behavior 
of the two-point correlation function of the fundamental fields. This asymptotic behavior is 
determined by the location of the pole of the two-point function, as in the continuum case. 
There is however a slight complication due to the absence of rotation invariance outside 
the scaling region. For sake of completeness, let us discuss this point, that was solved in 
the 0(A) case by Miiller, Raddatz and Ruhl [ |139| ] by applying saddle-point techniques to 
estimate the large- distance behavior of G(x) (see also Cristofano et al. |140|| ). 



For sake of definiteness and simplicity we shall now focus on the case C2 = 0, where the 
inverse two-point function is 



1 



[11.1) 



Let us introduce the function fJ,(9), indicating the coefficient of the large-distance exponential 
decay, which is dependent on the polar angle 9 in the (x, y) plane, and define 



yu(6>) = fi x (9) cos 9 + fiy (9) sin 9, 
where \x x and \i v are functions of 9 such that 

G' 1 (pi=ifii) = 0. 



;ii.2) 



;n.3) 



Eq. ( p.l.3| ) is not sufficient to determine /J,i(6), but, guided by the saddle-point analysis of 
Ref. |p.39|| ) we make the large- A Ansatz 

(11.4a) 
(11.4b) 

(11.5) 

p 2 + m 2 Q = (11.6) 



lMi = tM* + ^l*ii + o(Jp 
sinh fj, 0x = u cos 9, sinh fi 0y = u sin 9. 
The generalization to C2 7^ would consist in the replacement 

sinh fj, 0i — > Ci sinh fi 0i + \c2 sinh 2fi 0i . 
We can now solve the equation 



for pi = ifiQi with the Ansatz (|11.4|) and obtain 



u = m 



2 + X 



1 + 



m 



cos 29 
A + ml 



-1/2 



;n.7) 



hence 
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yU (#) = cos # arsinh(z/ cos #) + sm axsinh^o sin 0) = m + 0(ml) . (H-8) 
Now replacing Eq. ( 11.4 a) in the condition ( 11.3Q we immediately get 

- 2 [uq/jlix cos 9 + z/ /% sin 0] + A(z// i) = 0, (11.9) 



implying 



m{0) = cos e + fi ly sin e = (n.i ) 



and 



/i 2 W = A*S(*) + ^ ^ + O (J^j . (11.11) 

These results are quite general and may be useful in detailed studies of lattice models 
outside the scaling region. One may check that in the appropriate (scaling) limit the de- 
pendence on 9 disappears in /ii(6 l ) as well as in fiQ. However, in practice we may limit our 
attention to only two special values of 9, corresponding to the two extrema 9 = (side 
correlation) and 9 = j7r (diagonal correlation). In both cases no special Ansatz is needed in 
order to find the expression for the mass-gap: in the first case the lattice symmetry implies 
fi y (0) = and therefore /i(0) = f-i x (0); in the second case we have fi x (^7r) = H y {\^) = \ia 
and therefore fJ>(jir) = V% fJ>d- Eq. (|11.3|) reduces in both cases to a single- variable equation. 
It is easy to check that, when C2 = 0, 

/i (0) = 2arsinhim , (11.12a) 

// ( |tt) = 2V2 arsinh ; (11.12b) 

2y 2 

/io(0) and /^q(t7t) are plotted as functions of / in Fig. [IB]. It is also pleasant to check that, in 
these two specific cases, the resulting definition of the mass gap coincides with that obtained 
from the so-called wall-wall correlations 



G s (y-x) = ^-^2,G{xx,x;y 1 ,y), (11.13a) 

xi,yi 

Gd = ~r ^2 G ( xi > x ~ Xi >y i >y~yd> (n.i3b) 



xi, yi 



whose effective one-dimensional propagators in momentum space are 

G- x = G- x {p x =Q,p v ), (H.14a) 
G, 1 = ■±=G- 1 (p x =p d ,p y = Pd ). (11.14b) 

Since side and diagonal correlations are the two extremal cases on the lattice, rotation 
invariance may be verified simply by checking that 
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G s (x) = G d {x) 



(11.15) 



for all values of x. 

As long as we are interested only in the scaling behavior, in order to minimize the effort 
we can focus on the case 8 = 0. We can be fairly general, because in the scaling region 

//o(0) =m + O{m 3 ), (11.16a) 
v (0) =m + O(m 3 ), (11.16b) 

independently of the choice of the lattice action. Therefore 

= ™o + ^ Ei(irno) + 0(mj) + O (J-\ , (11.17) 
and we need only to compute the lattice 0(1/ N) contribution to the self-energy. 

B. Scaling behavior of the free energy 

Before attacking the above-mentioned problem, let us however notice that, if we believe 
standard renormalization-group arguments, we may extract the scaling behavior on the 
lattice directly and with lesser effort from the expression of the lattice free energy. The 
subtle point in the evaluation of the lattice free energy has to do with the existence of 
a perturbative tail that obscures the scaling behavior of this physical quantity. However 
this turns out to be a tractable problem, because the asymptotic expansion introduced in 
Sect. H allows us to isolate and remove unambiguously the perturbative tail of the free 
energy. In practice, without belaboring in this point (cf. Ref. |J for more details), we must 
simply translate Eq. ( |5.2[) into the lattice language and write, in analogy with the continuum 
notation, 

F (scal) = N Tr ln J -2 + m 2) _ N Tr ln -2 _ Nf3m 2 

+ iTrln^^-iTrln^^L^ 

+ ^^^;^(p)-|TrlnZ\^ M »| m2o=0 , (11.18) 
where we need to recall that 

= ^ + A le]iP) C(P) ( n - 19 ) 

and 

THnZl-i (p) = / lndetZU (p). (11.20) 
Hence we only need to evaluate 



11.21) 



52 



The mo — > limits of the relevant expressions are obtained by the replacements 



(9) 
) 



C 0li + PD, 



0/i ■ 



(11.22a) 
(11.22b) 
(11.22c) 



Let us now notice that, according to the rules of the asymptotic expansions, we may 
evaluate 



N Tr In (p 2 + m 2 ) - N Tr In p 2 - N(3m 2 
d 2 p ml AT/D 2 , , r f d 2 p 



N 



(2tt) 2 p 



^ - Npm 2 + N 



(2tt) 



In 1 + ^ 



Ulr, 



P 



0(m 4 Q ) 



^-ml + 0{mt). 



11.23) 



Moreover, if we treat TOq and (3 as formally independent variables in Eqs. ( p.4|) and ( |9.33|) , 
we can introduce partial derivatives via the relationship 



dm 2 . 



d dp d 







dml dm 2 , 9/3 



dm 2 ) 



0/? 



11.24) 



this allows us to reformulate Eq. (|11.18|) , obtaining the following representation of the scaling 
lattice free energy to 0(1/N): 



F (scal) = il m 



N 

47T 



d 



- m, 



d 2 p 
(2tt) 2 



d 2 

-^JLlndetA-)^) 



mg=0 



continuum counterterms. 



11.25) 



The continuum conterterms in Eq. ( |11.25|) are to be introduced according to the usual rules 
and will generate a contribution akin to the counterterms used in the sharp-momentum cutoff 
evaluation of the free energy. We remind that the continuum conterparts of the relevant 
quantities appearing in Eq. flll.25| ) are 



i 



l n i±l + (a 2 ) 



2np i i i-\ 

2tt £ - 1 



+ 0(a 2 



l_ f d 2 q c 1 g 2 + 4c 2 E,g'cos 2 ^ 
4/3 7 (2tt) 2 c^ + c^ql^hv 



Co 



(11.26a) 
(11.26b) 
(11.26c) 



C - ComI p =o- 



The form of the continuum counterterms is therefore 



d 2 p 



In 



1 



2^-^-1 



In 



In 



2tt £-1 



11.27) 
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We confirm the observation that, in order to get a better agreement with the corresponding 
continuum result, it is convenient to redefine x by 



x 



X 



11.28) 



where k is to be kept constant along the renormalization-group trajectory. The redefinition 
( |11.28j ) is not going to affect Eq. ( |11.25[ ), because x turns out to be a perturbative function of 
j3 and is therefore not affected by the partial derivative with respect to TOq. Nonperturbative 
redefinitions of x are not allowed; as we shall see, they would spoil the scaling properties of 
the physical quantities. 

Eq. (|11. 2 71) can be analyzed in the light of the results presented in Sect. |, and we 
recognize that it can be rephrased in the form 



implying 



— ^ [ln4vr/? + (3 - 2vrx) ln(4vr/5 + 2vrx - 2) + c F {k)\ 
An 



An p 2 



+ 



3 - 2ttx 



ln(j> 2 / m o) ln(p 2 /ml) + Inn — 2 



_p(scal) _ _p(SM) + §p 



11.29) 



(11.30) 



where 



6F 







2 m ° dm, 



o 



L dp 2 An 
Air p 2 



d 2 p 

WY 2 



tJ^2 lndet4 (i» 



mg=0 



3 — 27TX 



ln(p 2 /m 2 ) ln(j9 2 / m o) + Inn — 2 
and F^ SM ^ is given by Eq. (|5.8|) (with x replaced by x) . 



11.31) 



C. Scaling behavior of the self-energy 



The 0(1/N) contribution to the self-energy is obtained by computing the Feynman 
diagrams drawn in Fig. pi. Evaluating the diagrams, according to the rules presented in 
Fig. O, leads to 



A(p) 



d 2 k Ar a) (k) 1 



+ 



(2tt) 2 p+k 2 + m 2 



[IV 



exp [\i{k^ - k v )] 



V^{p+\k,k)V Zv {p+\k,k) 

r 2 9 

p+k + TOq 
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+ 2A«)(°) 



d 2 k 

(2vr) 2 



^Z\ (e)Atl ,(fc) 



A" 7 



X 



dm 2 . 



^» u {k)-2{l + xf)5» v J 



d 2 g V^(g,A;) 



(2vr) 2 g 2 + m 2 

+ ( 1 + x/) A «,(0) (11.32, 



//// 



Let us now apply the decomposition of the total derivative with respect to the mass into 
partial derivatives, perform some simplifications, and rearrange the terms into the form 



d 2 A; 



A«)(*0 



(2tt) 2 

+ ^ (Q) (o) y 
+ 



1 d 



d 2 A; 



(2tt) 2 dml ^ 

3 ^A*)A»( fc ){ V^i(P>*0 



(27T) 

- exp [§i(/fy - fc ; 
d 2 A; 



+ ^ (Q) (o) y 

+ x/ 



(27T) 



^(?+§M)^(p+|M) 1 d A-i (k)6 ^ (k) } 



r z 9 

P+K + Wig 



2 



9m 2 lndet^^CAr) - 5>W/«<(*) ^^W^a 



d 2 A; 



(2tt)' 



(11.33) 



Eqs. ( |7.10| ) imply the following relationship: 



2^ exp (-2^ M j S ^ V^W, «) - exp ^HS - MJ ~ ; 



exp (—§«&„) 



7)2 _i_ 77)2 1 

° V 3v {p+\k,k)--(V 3v {p+\k,k) + V 3v {p-\k,k)) 



p+k + T71q 



(11.34) 

Eliminating terms vanishing under symmetric /c-integration, we can perform the replacement 

V^{p+\k,k)V 3v {p+\k,k) 



V^VGo, k) - exp [-z^ - fej 



W(p, *) SL (k) + 



i , , , P 2 + '"o 



p+/e + rriQ 



V 3 „{p+\k,k) 



V 3p (p+\k,k) 



(11.35) 



where 
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W( P ,k) = J2v t Ap,k)-j: lV — «**'* )] ' - (11-36) 
M M P+k 

In order to evaluate the function Ui(p) in the scaling region, we must keep in mind that 
the formalism introduced in Sect. |9] allows for a direct expansion of the lattice integrals in 
terms of the "external" variables m and p, since all the nonanaliticity is accumulated in the 
continuum integrals, which in turn have already been computed in the sharp-momentum 
regularization scheme. We now realize that, as far as the lattice contribution is concerned, 
we can make explicit use of the following relationship, holding for any regular (analytic) 
function /: 

f(k; p, m ) = f(k; 0, 0) + p 2 ^f(k; 0, 0) + m\ ^f(k] 0, 0). (11.37) 

In particular, when imposing the on-shell condition p 2 + m\ = 0, we find 

/(*; tm 0l mo) = f(k; 0, 0) - m 2 ^f(k; 0, 0) + m 2 ^f(k] 0, 0), (11.38) 

and, as a consequence, 

f(k; p, mo) = f(k; tm , m ) + (p 2 + m 2 ) -^f(k; 0, 0). (11.39) 

Applying Eq. ( |11.39| ) to Si(p), we come to the conclusion that 

£i(p, m ) = E[ s M) + 5m\ - (p 2 + m 2 ) 5Z 1 , (11.40) 

where X'j SM ' ) is the value of the self-energy in the sharp-momentum cutoff regularization 
scheme, and 5m 2 , 8Z\ are the 0(1/ N) lattice contributions to mass and wavefunction renor- 
malization, and are amenable to finite lattice integrals whose infrared regularization is pro- 
vided by the appropriate continuum counterterms. 



D. Lattice contribution to mass renormalization 



The 0(1/ N) contribution to the mass, (for arbitrary values of the coupling constant), 
can be obtained from Eqs. flll.llj ) and ( |11.33| ): 



/4 = — Ei(ifio) 



/A) 



d 2 k 



A {a) (k) 



1 



+ 



d 2 k 
(2vr) 2 



J2Ae),Ak)SlM 



k+ifl Q + !71q 

W(i^k) + l -A {a) (Q)^- 2 A^(k) 



+ Xf I ^EA^WV V^o,k)+A (a) (0)^^A^ )fM (k)) 
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^0 



+ 



d 2 k 
d 2 k 



1 9 



fc+^ 2 + m 2 2 



W{ino t k)-~A^(k) 



/d 2 k T (9 / 



d 2 A; 9 



(27r) 2 <9m 
d 2 £; 



2 [iB^C^+lBdet^ 1 ^) 

<9x 



(2vr) 2 



dm 2 . 



11.41^ 



In order to evaluate £m 2 , let us first complete the lattice computation of m\ in the 
scaling region from Eq. Q11.17 ). From the definition of the asymptotic expansion presented 
in Sect. |9] it is easy to check that, under symmetric A;- integration, the following replacements 
are allowed: 



1 



k+ifiQ + tuq 



B<?\k) + mlB<?\k)] + 0{ml) = \ ^~ a )(k) + 0(rn*), 



1 

dp 

We make now the crucial assumption 



W(ifio, k) 
V^(ifi , k) 



l^A 7e ](k) + 0(mt), 



0. 



dx 



(11.42a) 
;il.42b) 

(11.42c) 



11.43) 



which corresponds to our choice of renormalizat ion-group trajectories and defines the class 
of theories we are going to generate in the continuum limit. The renormalization-group flow 
trajectories defined by Eq. (|11.28|) , consistently with Eq. ( |11.43|) , are plotted in Fig. [L7| (for 
c 2 = 0). Under such assumption and making use of the identity 



we can conclude that 



Ex{imo) = 4:11ml x 



4 (a )(0) =47rm 2 + OK), 
d 2 A; d 



11.44) 



(2tt) 2 dm 2 



ln^AO+lndet^^A;) 



mg=0 



+ O(m ) + continuum counterterms. 



11.45) 



In order to write down explicitly the continuum counterterms appropriate to Ei{imo), it is 
convenient to go back to Eq. ( |11.32j ) and perform the replacements indicated in Eqs. ( |11.26|) , 
supplemented by the substitutions 
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1 : z 9 

k+vm + tjIq 
W(im , k) 
V 4tl (im , k) 



1 

i. 



The resulting expression for the continuum conterterms is then 



(11.46a) 

;il.46b) 
(11.46c) 
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1 - T 



(mo) 



2tt 



111 



2tt 



£ln4^ - 2 + 2nx 



£-1 



11.47) 



Let us notice that Eq. (|11.47 ) was derived from Eq. flll.32 ) without making any assumption 
on the dependence of x on (3. 

In the light of the results of Sect. |5], Eq. ( |11.47| ) turns into 



m5 
Air 



[\nAn(3 + (3 - 2ttx) ln(47r/3 + 2vrx - 2) + c m (x) 



M Z dk 2 4tt 



4tt k 2 



+ 



3 — 27TX 



ln{k 2 /m 2 ) \n{k 2 /m 2 ) + 2vrx - 2 



implying 



where 



2(SM) 



+ 6ml, 



Sml = Air 5F 



(11.48) 



11.49) 



11.50) 



and m!f M ^ is given by Eq. (]5.22| a) (with x replaced by x). 

In conclusion, we can summarize the results contained in Eqs. ( 11.30 ) and ( 11.50 ) into 
the relationship 



= [N- c m (x) + CF (x)] ^ + O (J^j + 0( 



m ). 



11.51) 



Eq. ( p.1.511) is one of the main results of the present Section. It shows that, performing a 
proper perturbative redefinition of the parameter x, it is possible to map the lattice theory in 
the scaling region into the corresponding continuum theory in such a way that the fundamen- 
tal scaling relationships are preserved. Eq. ( |11.51|) is identical in form to the corresponding 
continuum relationship ( |5.29| ). Let us notice once more that nonperturbative redefinitions 
of x would spoil the scaling relationship, because they would change F^ scal ) without affect- 
ing Si (the partial derivative dx/dm 2 , cancels in Eq. ( 11.41| ) with the corresponding term 
coming from d ^ v / dmfy . 
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Recalling the explicit form of the asymptotic expansions presented in Sect. |9], we can 
now construct the following representation of the 0(1/N) contribution to the physical mass 
gap in the scaling region: 

d 2 k A[ a) +(3B[ a) 



5m 1 = litm, 




4 Q) + PB$ 



(o) 



+ 



(Coi + /?A)i)(Ci 2 + /?£>i 2 ) + (C 02 + PD m )(C n + /3Du 



X 



x 



E(Co,+^)|(4 e) +^r 

. A 4 



C + /3 



idfc 2 2 



47r A; 2 



+ 



3 — 27TX 



ln(P/m2) ln(P/m2) + 2vrx-2 



11.52) 



The integral representation is improper, since the denominators vanish on certain curves; 
the generalized principal-part prescription giving a precise meaning to Eq. ( |11.52| ) will be 
discussed in Subsect. |11 G| (cf. Ref. ||). 

A rather straightforward consequence of Eqs. ( p. 1.37 ) and ( |11.38| ) is the relationship 



d 

f(k; im , m ) = f(k; 0, m ) - m 2 T^f(k; 0, m ) + 0{m 



1 
0)1 



11.53) 



implying that the lattice contribution to Zi(imo) and to ^i(O) — m^Z^O) have the same 
scaling limit. As a consequence we should have obtained the same expression for 5m 2 if we 
had computed the lattice counterpart of m\ in the scaling region. 

We are now ready to apply our knowledge of Sm 2 to the problem of finding an explicit 
expression for any physical quantity with known scaling properties, when we analyze it on 
the lattice in the scaling regime. Let us assume Q to be a physical quantity with canonical 
scaling dimension 2Sq in mass units (i.e. no anomalous dimension), and let us compute its 
large N limit and 0(1/N) corrections: 



Q = (miy 



N 



N 2 



11.54) 



evaluating the 0{1/N) corrections in the SM regularization scheme. Since go is universal 
and /^-independent, we can immediately predict the lattice scaling behavior of Q to be 

1 M , ^( I 2\<5q+2 



Q 



K) 



2\ 5 Q 



1 ... 1 . Sm 2 (P) 



TTlr 



+ 1 (ml) 



11.55) 



Eq. (|11.55|) realizes the promise of evaluating every physical quantity in the scaling regime 
by means of a single lattice computation. 
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E. Lattice contribution to wavefunction renormalization 



The discussion of the wavefunction renormalization goes along the same lines. We recall 
from Eq. ( |L3g ) and fllLlOD that 



5Z 1 







dp' 



p=0 
mo=0 



11.56) 



From Eq. ( |11.33| ) we obtain 
-8Z X = 



d 2 k A d 1 

A (a) {k) 



(2vr) 2 
d 2 k 



d P 2 p+k 



p=0 



J ^ A m^k)8lM^w{p,k) 



p=0 



(2vr) 2 
d 2 k 



(2vr) 2 

where we made use of the relationships 



A(e)^ u (k) 8^ u Di^(k) + continuum conterterms, 



(11.57) 



V^k) = V ^f' k) =Do„(fc), 



d 



p=0 



-D^{k) . 



(11.58a) 
;il.58b) 



Let us further notice that 
d 1 



dp 2 p+k' 
d 

■W(p,k) 



dp- 



p=0 

P =o 2 



(k 2 

lB[ e \k)-^j2K{ D o^)] 2 . 



(11.59a) 
(11.59b) 



We are finally ready to write down the complete explicit expression for bZ\\ 
d 2 fc 1 



SZ X 



(2tt) 2 A^ + PB { a 



(k 2 Y 



+ 



d 2 A; 



k 2 



k 2 



i 



k; 



(C 01 + /3A,i)Am + (C 02 + /?A) 2 )Ah - V(C 0At + PD 0ll )-£D 0v 
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+ ^ ~ ~~^7" [(0)1 + /?A)l)£>12 + (G)2 + ^02)^11] 



2 p+a 







+ 



f £ dA; 2 



47T 



2tt 



47T 



1 

+ - 



In (k 2 /m 2 ) In (fc 2 /m 2 ) + 2vrx - 2 k\ k 2 



11.60) 



and we may verify that 



lim 6Zi = 0. 



11.61) 



F. CP^ -1 and 0(2iV) models 

At the borders of the parameter space, CP^ -1 (f3 v = 0) and 0(2N) (/3 g = 0) models 
require a separate discussion, because of some subtleties related to the order of the limiting 
procedures in the lattice formulation. 

We know that, because of gauge invariance and confinement, 



lim mi = 00. 

i*->0 



(11.62) 



However, the quantity 5m\ stays finite, and this is important because of its renormalizat ion- 
group interpretation: it allows us to define a renormalization-group-invariant, l/N- 
expandable scale in lattice CP^" 1 models. When computing 5m\ for CP^ -1 models, it 
is necessary to keep in mind the transversality of the vector propagator, due to gauge invari- 
ance. As a consequence, one may extract directly from Eq. ( |11.41| ) the correct expression 



5m iU=o 



2nm 2 ] 



d 2 fc 



A{ a) + (3B\ a > A*> + (5B\ 
A ( Q a) + (3B ( Q a) + 



?(«) 



(0) 



A 



dk 2 2 



47T k 2 



+ 



ln(fc 2 /mQ) ln(fc 2 /m, 



11.63) 



Eq. Q1T63D differs from the k -> limit of Eq. ( [jXggj ): 



d 2 fc 
(2V0 2 



?(<*) 



+ 



l(«) 



+ 



^dA; 2 2 



47T /C 2 



ln(/c 2 /mQ) ln(fc 2 /mo) 



11.64) 



On the other side, when we consider 0(2iV) models we find from Eq. Qll.41 
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and 



5Z X 



/3 g =0 



d 2 k 



(2vr) 2 A^+(3B ( 



(a) 



(a) 



M 



Zdk 2 1 



2tt 



4tt k 2 In (k 2 /ml) ' 



;il.65b) 



Eqs. ( ll.65| ) ought to be compared with the x — > oo limit of Eqs. (11.52) and ( |11.6U ): 



lim 5m\ 



271ml 



+ 



d 2 k A^+fiBl 



(a) 



£ dfc 2 2 



(2tt) 2 A («) + p B to J 
(2vr) 2 ^C 0ll + f3D J 



47T fc 2 



In (k 2 /ml) 



11.66a) 



and 



lim SZi 



+ 



d 2 A; 



2tt) 2 
d 2 A; 



.4, 



(«) 
o 



2 1 



(a) 



+ 



(k 2 



E 



tdk 2 1 



2tt 



47r P ln(fc 2 /mo) 



;il.66b) 



We may observe that Eqs. ( |11.63| ) and ( 11.65|) differ from Eqs. ( |11.64| ) and ( |11.66| ) respec- 
tively because of a term depending on C*i M and Di^. Therefore this difference would be 
completely absent in the SM and in all other continuum regularization schemes. On the 
lattice these contributions are nonvanishing; however they are perturbative in 1//3, and in 
particular when k — > oo they do not affect the scaling relationship ( |11.51| ), where 



lim [c m (x) - cp(x)] = 1. 

rc— >oo 



11.67) 



Therefore in the scaling region they are amenable to a perturbative redefinition of the 
coupling /, i.e. to a different regularization scheme of the same theory. We shall come back 
to this phenomenon in Sect. 14, where we shall explicitly discuss the (3 — ► oo limit of the 
above expressions. 



G. Evaluation of lattice integrals 

The numerical evaluation of Eqs. ( |11.52| ) and (|11.60| ) is a fairly nontrivial task. The 
general structure of Sm 2 and SZi is that of a difference between an infrared- singular lattice 
integral and an ultraviolet-cutoffed continuum integral, with the same singular infrared 
behavior as the lattice integral. 

By replacing ln(k 2 /ml) with An (3 + ln(/c 2 /M|), it is possible to perform the expansion of 
bm\ in a power series in the powers of / = 1/(2/5). The coefficients of this weak-coupling 
expansion are individually infrared-regular combinations of lattice and continuum integrals. 
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It is easy to generate the numerical values of the coefficients up to high orders, simply by 
expanding the integrands and integrating term by term. However the full expressions them- 
selves are not proper integrals, reflecting the fact that the series are not Borel-summable. 
As a consequence, the approximation by series expansion will fail for sufficiently large values 
of/- 

In practice, we found that numerical predictions from series evaluation become unstable 
(with respect to different truncations) at the border of the scaling region. We don't mean 
that perturbative evaluations be intrinsically meaningless; we only stress that extrapolation 
to the intermediate coupling may be dangerous and need to be carefully tested against 
instabilities. In particular, disagreement between Monte Carlo results and low orders of 
perturbation theory cannot be automatically interpreted as absence of scaling and/or failure 
of field theoretical predictions. 

In any case, in the models we are discussing an explicit way out of these difficulties can be 
found, since a generalized principal-part prescription allows an unambiguous evaluation and 
resummation of (the scaling part of) the series (cf. Ref. 0). Our prescription is expressed 
as follows: 

1. represent the integrals as sums of individual terms of the form 



d 2 k 



M(k) 



(2vr) 2 [(3 + N(k)] r 



2. define complex variable functions 



3. identify 



d 2 k 



C{z) 



M(k) 



dz' 



z — z 



d 2 k 



M(k)S(z' + N(k)); 



:-v, 



(2vr) 2 [(3 + N(k)f 



(n- 1)! d/3 



n-l 



lim 



C(/3 + ie) + C(P - is) 



11.68) 



11.69) 



11.70) 



From a numerical point of view, a direct principal-part integration is possible but very 
unstable, since the line of vanishing denominator of the integrands has to be found numer- 
ically, and moreover it shrinks rapidly towards zero in the limit / —>■ 0. The solution is to 
adopt the method described in Ref. 0: the integral is split into the two regions k 2 < p 2 
and k 2 > p 2 , choosing p such that all the singularities are included into the first region. The 
integral over the second region is perfectly regular and it is easily computed numerically. 
In the first region, we expand the integrand in a power series in k and integrate analyti- 
cally term by term, resulting in the exponential-integral and related functions. The typical 
analytic integration involved in the abovementioned procedure has the form 



Jo 



dx 



x" 



(hix 



8+1 



(n + 1) 



si 



y n+1 S s ((n+l)lny) 



11.71) 



where 



t\ 

yt+l 



11.72) 



t=o 



t=s 
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and Ei is the standard exponential-integral function. Tuning p as a function of /, the 
nth order of this expansion approaches the exact result with an error decreasing as fast as 
exp[— 4ir(n + l)/3\. The series expansion is shown numerically to approximate extremely well 
the integrand in the relevant region. 

It is worth noticing that the resummation procedure we have sketched turns out to be 
useful also in the evaluation of the so-called perturbative tails of composite operators, i.e. 
those non-scaling contributions to vacuum expectation values resulting from spin waves that 
must be explicitly subtracted in order to derive normal-ordered quantum expectation values 

§• 

We may notice that the ambiguity in lim^^ C(z) is related to the Borel ambiguity in the 
resummation of the perturbative series observed by David fl80U8l| . However the coefficients 
of the asymptotic expansion in powers of wig of a physical 1/iV-expandable quantity like 
the mass gap are real functions, whose expansion in powers of the coupling is the standard 
perturbative series. 



12. EVALUATION OF PHYSICAL QUANTITIES FOR NEAREST-NEIGHBOR 

INTERACTIONS. 



As already discussed at the end of Sect. ||, the task of performing accurate numerical eval- 
uations of the integrals entering the 0(1/N) contributions to physical quantities is greatly 
simplified in the case of nearest-neighbor interactions C2 = 0, where we possess integral 
representations of the lattice propagators. 

Notable simplifications occur in Eq. ( |11.41|) in this special case: the exact 0(1/N) con- 
tribution to the mass is 



A*? 



^0 



+ 



d 2 k 

d 2 k 

(2*0 



a A«)(*) 



k+ifj,o + m 



^(a)(*) 



1 Aa)(0) 



A 1 



o( 4 + m o)- 



k+2ifi 



k+ip, + ml 



+ 2A«)(°) 



dmg 




where 



.1 1 mi 
A' 1 = 1 1 

4(3 + 4 ' 



x 



x ■ 



4(3 



4(3-1 



(12.1) 
(12.2a) 

12.2b) 



and A( a \, Arg) can be computed from Eqs. ( fj,19| ). 

Eq. ( |12.1|) can be tested, for sufficiently small /3, against existing strong-coupling results. 
Strong-coupling expansions for the mass gap in 0(N) and CP^^ 1 models were pioneered in 



Refs. ||131| , |132| , |141| , |142|| , and in particular it was recognized that, in the CP case, strong- 



coupling and large- N limits cannot trivially commute: indeed the U(l) gauge invariance is 



64 



spontaneously broken at N = oo (states in the fundamental representation are deconfmed 
free particles); on the other hand, it cannot be broken in strong coupling for any finite 
value of N. A signal for this phenomenon is the nonexistence of a regular strong-coupling 
expansion of Eq. (|12.1|) at x = 0. In the 0(N) case, extended strong-coupling series for the 
mass gap (up to 0(/3 n )) were obtained by Butera and coworkers in Ref. [|143|| and analyzed 
in the large- N limit in Ref. |[L44j| . Results up to 0(/3 14 ) for the magnetic susceptibility 



and related quantities were obtained in Refs. [ |145] , |146| | . Strong-coupling expansion of O(N) 



models with improved action was discussed in Ref. [147 



Eq. ( |12.1 ) can be further simplified in the large- x regime, where one can perform a 
1/x expansion of the integral involving the vector field propagator, leading to an expression 
which can be evaluated exactly to 0(l/x). We obtain the result 



Mi 



d 2 k 



(27T 



„A«)(*) 



i 



x 2 1 4/3 



2 2 
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(0) 



(a) 



(k) 



dm 2 . 



+ 



1 A«)(°) 



A- 1 



112.3) 



Eq. ( |12.3| ) is a rather good approximation of Eq. ( 12.1 ) even for quite small values of x, as 
long as we consider not too large values of /3. Eq. ( |12.3| ) should be compared with the exact 
0(1/ N) contribution to the mass of the 0(2iV) model 




A {a) {k) 



' + 2^)(°)^ 



k+i/io +m„ 



;i2.4) 



When C2 = significant simplifications occur also in the scaling region contributions 
obtained in Eqs. ( |11.52|) and ( |1 1.60| ) . The explicit forms of 8m\ and 8Z\ are respectively 



5m\ = 27rm, 




d 2 k 
32 dP 2 



+ (3B\ a> A { °> + /3B\° ' + /3x/(4/3 - 1) 
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4 e) + /3£^ + x 

3 — 27TX 



47T fc 2 
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4/3 



ln(fc 2 /m 2 ) ln(fc 2 /m 2 ) + 27rx-2 
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and 



8Z X 



d 2 k \Br + \/{p 



(2vr) 2 A^+PB^ 



1 



2 4/3 



+ 



1q -r p-l^ 

" 32 dfc 2 : 



d 2 A; B[ e) + x/(4/3 - 1) 
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f 2 J (2vr) 2 A® + p B <P + Z 
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2tt 



47T 



4tt A; 2 Lln(fc 2 /m 2 ) In (P/m 2 ) + 2vrx - 2_ 



;i2.6) 



It is easy to compute the first nontrivial orders in the weak coupling expansion of Eqs. (|12.5[) 

2) + (c! Q) + cf + x) / + O (f 2 ) , 



and ( |12.6| ), obtaining 

5m\ 



3tt 

T 



12.7) 
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and 



f 

8Z X = — 
2tt 



3tt 



^-2) + (cr+cr + x+±) / + 0( / 2 ) 



7T 



where 
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7T 
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v4q 5 



5, 



d 2 A; 



(a) 




(a) 



5^ 5, 







2ttP 



5, 



jgW fl W + 16 



2ttA; 2 



We can now make use of a number of identities holding for regulated lattice 
discussed in Appendix [F], to obtain 



8 2tt 1 



7T 



8 2tt 1 



0.0282552, 



1.7120726 



(12.8) 

(12.9a) 

(12.9b) 
integrals, 

(12.10a) 
(12.10b) 



(cf. Eq. (|F3p). 

Eqs. ( |12.7 ) and ( |12.8| ) have a straightforward relationship with the three-loop computa- 
tion of the lattice (3 and 7 functions. In particular the 0(2iV) results 



(12.11a) 
;i2.11b) 



reproduce the original calculation by Falcioni and Treves [148], already confirmed in 



Refs. |TT6| , |149 1. 

Higher-order coefficients of the weak coupling series can be evaluated numerically with 
high precision, and different truncations of the series can be compared, checking for stability. 

We can also compute numerically to high precision the difference between the exact 
lattice representation ( 12.1|) of /i 2 and its SM continuum counterpart m\ SM as expressed by 
Eq. ( |5.24| ) , in the region /3 < 1.5. 

Finally we may evaluate the representation ( 12.5Q of the scaling contribution to the 
mass gap Smf, as well as the representation (12.6) of the renormalization function 8Z\. We 
evaluated the integrals in Eqs. ( |12.5|) and ( 12.6 ) using the expansion in exponential-integral 
functions described in Sect. |Tl| and Ref. p. 

These different evaluations should all agree with each other in the very weak-coupling 
domain, where truncated perturbative series are accurate. Moreover the difference between 
the last two determinations of 8m\ is entirely due to scaling violations; therefore it can 
be compared with independent determinations of the scaling region, such as the study of 
rotation invariance properties of the mass gap. 

All the relevant numerical results are presented in Figs. [18] and |19, where 8m\ and 8Z\ 
respectively are plotted as functions of / for different values of x. 
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In the case c<i 
simply 

5m\ = 2-Kml 



we can also perform a 1/x expansion of Eq. (|12.5|) . The result is 



d 2 k A[ a) +f3B[ 
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(2tt) 2 A^+f3B t 
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2(3 1_ 
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7T 4 4/3-1 
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4tt 2 



X 



12.12) 



It is possible to check directly the consistency of Eq. ( |12.12|) with the asymptotic expansion 
of Eq. ( |12.3| ). As a byproduct, one may also verify that the asymptotic behavior of c m (x) 
is correctly represented by Eq. Q5.26]) . The corresponding result for the 0(2N) models is 
expressed by Eq. ( |11.65| a) and is plotted in Fig. |2D|. 

In the opposite limit, the result for CP^ -1 models as expressed by Eq. ( 11.63 ) is plotted 
in Fig. EH. 



13. TOPOLOGICAL OPERATORS 



We already mentioned at the end of Sect. |5] the special role played by topological prop- 
erties in the limit x = 0, corresponding to pure CP^ 1 models and U(l) gauge mvanance. 
The problem of defining a sensible lattice counterpart of the topological charge density (|5.66|) 
has long been debated in the literature. 



The geometrical definition originally proposed by Berg and Liischer [ |150|| amounts to 
defining 

q& = — Im{lntr(P n ) + \ntr(P n+u P n+fJi+u P n )} , p^i/, (13-1) 

Z7T 

where 

Pn,ij Zn,iZn,j • (13-2) 

g| has the advantage of generating integer values of the topological charge for any given field 
configuration, and one can prove the absence of a perturbative tail in the 1/N expansion 
134| |. In formulations involving an explicit U(l) gauge field X n ^ = exp^ n /i , an alternative 



geometrical definition is obtained by defining 

^ ^ £fiv(@n,[i "I" @n+Li,v @n+u,fM @n,v) ■ (13.3) 



47T 



q n enjoys the same properties of g§ and has the same 1/N expansion. 

Unfortunately for finite N geometrical definitions are plagued by the so-called "disloca- 
tions", and therefore one cannot extract the correct scaling behavior from numerical data 



151| , |152|| . It is however possible to express the topological charge in terms of a local operator 
constructed from the gauge fields. Let us define the plaquette operator (elementary Wilson 
loop): 
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u 



«n = «n,i2 = u n,2i = exp{27ug n }. 



(13.4a) 
(13.4b) 



Taking proper combinations of higher powers of the plaquette operator, it is possible to 
construct an infinite sequence of local operators 



Im {(«„)'} , 



;i3.5) 



whose formal k — > oo limit is exactly Eq. ( [13. 3j ). One can then construct the sequence of 
topological susceptibilities 



Xt 



;i3.6) 



The perturbative evaluation of these quantities is obtained by considering that q n is 
linear in the effective Lagrangian field #„ jA1 , and expanding in a power series in q n . One 
may easily show that 



2vrgf = 2irq n - 



(27rg n ) 2fc+1 + 0(£+*) 



2k+3\ 



(2fc+l)l 

In standard perturbation theory, using lowest-order momentum space propagators 

(g(p)g(-p)) = ^p^p 2 , 



;i3.7) 



;i3.8) 



one can prove the relationship 

(q {k \p)q (k) (-p)) = (q(p)q(-p)) 



sowar 



A more refined analysis, based on the relationship 

X[ k) = lim (q^ k \p)q {k \-p)) 



and the observation that 



lim (q(p)q(-p)) = 0, 



o 2 ^0 



;is.9) 



(13.10) 



(13.11) 



allows us to prove that, in the leading order, 



i / • i / 2/ 

Xt - Xt + c fc ( — 



2k+l 



13.12) 



where 
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(13.13) 



for large k. Eq. (|13.12[) shows that the perturbative tail of the topological susceptibility 



involves for high k only very high powers of /. The same property might easily be shown 
to hold also for other mixing coefficients. However, the corresponding numerical weights 
are growing so fast with k that the convergence to the geometrical definition Xt cannot 
be uniform, i.e. the limit k — > oo does not commute with the continuum limit (/ — > 0). 
This phenomenon leads to a perturbative explanation of the observed discrepancy between 
geometrical and local-operator definitions of the topological susceptibility, and shows that, 
for fixed f /N, an optimal value of k should exist such that the mixing is minimized. 

From the point of view of the 1/N expansion, the situation is however quite different: for 
fixed / the perturbative tail is depressed by a factor (1/N) 2k , and therefore the absence of a 
perturbative tail of the geometrical definition in the 1/N expansion is confirmed; the differ- 
ence xi ~ Xt is calculable order by order in 1/N by generalizing the techniques described 
in the previous Sections, while Xt itself is simply obtained from the lattice counterpart of 
Eq. (ET57D : 



X t = ^ o j^y 2 fA {e) ( P ), (13.14) 

where A^(p) is the full lattice propagator of the field 9^. As a matter of illustration, let 
us consider the first contribution to the difference Xt ~ Xt, drawn in Fig. 22. A factor 
of j3 2 Z\(0)(p) is associated with each wavy line; as a consequence, the infrared behavior is 
regular. 



14. RATIO OF A PARAMETERS AND RENORMALIZATION-GROUP 

FUNCTIONS 



It is possible to analyze the results presented in Sect. [TTJ, and especially Eq. ( |11.52|) 



from the point of view of the perturbative renormalization group. Let us focus on the 
contributions to mf that depend on the specific lattice model adopted, i.e. the quantity 8m\ 
defined by Eq. (|TL49| ). 

A. Ratio of A parameters 



Let us consider the (3 — > oo limit of Eq. (|11.52| ) and notice that even in this limit a 
(x-independent) contribution to 8m\ survives: 



lim 8m\ = 2nm^ I 

/3^oo J 



d 2 k 

W) 2 



(14.1) 



Because of the noncommutativity of the limits, we must consider separately the borders of 
the parameter space. Indeed at x = we find 
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When x — > oo we must face an even more involved situation, because the limits x — > oo and 
/? — > oo do not commute, as one may easily check directly from the standard perturbative 
expansion. From Eq. ( |11.66| a) we obtain 



lim lim Sm 2 

(3 — >oo >c — >oo 



27rmg 



d 2 k 
(2vr) 2 



r( q ) n 



B 



(a) 



Finally when k = oo we obtain from Eq. ( |11.65| a) 

d 2 k 



lim Sm-, 



I3 g =0 



27rm n 



(2k) 



t i 



B[ a) 2_ 

B (a) + 12 



2 



14.3) 



(14.4) 



We observe that the difference between different regularizations of the same physical models 
is amenable in this limit to the condition Di^ ^ 0. 

Let us now come to the physical interpretation of Eqs. (|l4l|) , fll4~2p , ([143D , and ( |i4l|) . 
These quantities are obviously related to the ratio of the so-called lattice A parameter to 
the continuum (SM) A parameter in the models at hand ||1 53| , |1 54|| . To be more precise, and 
recalling Eqs. ( |9.19| ) and ( |9.20| ), we may write the relationship 



A 



SM 



1/ 



Mr 1 



5m\ 

2N p^'oo ml 



lim 



14.5) 



Actually we can do better than Eq. (|14.5|) ; we can exploit the fact that the ratio of the A 
parameters is essentially a one-loop phenomenon and our knowledge of the first coefficient of 
the renormalization-group (3 function (N for all x ^ oo, N—l when x — > oo) to exponentiate 
Eq. ( |14.5[ ) and obtain the exact relationships 
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when x ^ and oo, and finally 

IsM 



A, 



Ml exp 



7T 



N 



d 2 k 



(a) 



B 



(14.8) 



when x — > oo. Eq. (|14.8|) is to be compared with 
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A 



SM 



A, 



M L exp 



7T 



N — 1 



d 2 k 



Bl 



(a) 



B, 



(a) + k 2 



14.9) 



which we would obtain by setting /3 g = from the very beginning (standard 0(2N) models). 

We can now obtain explicit representations of the quantities entering Eqs. ( |14.6|) , ( |14.7|) , 
( P-4.8Q , and (|14.9| ). First we notice that 



X]« Dip k 1 . 2 1t/' a 1 1 r \ 

^ ~ f 2 = T^}^ Sm 2*7* ( C l + 4C 2 COS £M > 



ci + 4c 2 cos i/c 



0/j 



4 ^ ci + c 2 cos 2 



(14.10a) 
(14.10b) 



Furthermore, using the results presented in Sect. we can show that 



B[ a) BJ' 



B, 
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B, 



(0) 
o 



B[ a) 2_ 



Ci + c 2 cos 2 i/c^ 



0a< 



(14.11a) 
(14.11b) 



We notice that terms proportional to d^d^ In /c 2 are total derivatives that can be inte- 
grated exactly for any physically acceptable form of k 2 : 



1 

7T 



(14.12) 



These terms are not really lattice artifacts: they are related to the ratio A SM /Aj^, where 
y^Hs is the A parameter defined in the dimensional regularization scheme with minimal sub- 
traction (notice that in dimensional regularization the integral of a total derivative vanishes 
exactly). Therefore we obtain 



A 



SM 



A 



exp 



MS 



A 



SM 



A 



exp 



MS 
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'iV 
1 



[x 7^ oo) 



2(iV 
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(14.13a) 
(14.13b) 



Eqs. ( |14.13| ) are crucial in finding the variable change from SM to MS scheme and verifying 
the perturbative consistency of the continuum results. The ratio of A parameters can now 
be expressed in the more conventional form 



^MS _ 

Al 

when h = 0. 



Ml exp 




2c x c 2 
k 2 



£ 



sin 4 -hp 
ci + c 2 cos 2 i/c M 



(14.14) 
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7T 
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Mi = Ml exp i - 7 



when x^O and x 7^ 00, 



2cic 2 
fc 2 



sin 4 



Ci + c 2 cos 2 ifcp 



+ 



0/i 



At. 



when x — > 00, and 



Ml exp 



N — 1 



d 2 k 
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, A- 2 
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(14.15) 



(14.16) 
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Ml exp 



7T 



d 2 k Y^D^kl 



N-lJ (2n)i E ^kl 



(14.17) 



in the pure 0(2N) case. 

When C2 = all integrals can be computed in closed form; we obtain 
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Ams 
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AmE 



A, 



32 exp 
32 exp 
32 exp 
32 exp 
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i(N-l) 



(CP^ 1 ), 
(x^O), 
(x — > 00), 
(0(2iV)). 



(14.18) 
(14.19) 
(14.20) 
(14.21) 



Eqs. (|14.14|) , ( P-4.15Q , ( |14.16| ), and ( |14.17| ) have been all explicitly verified in standard 
perturbation theory. In particular, the difference between Eq. (|14.14j) and Eq. (|14.15|) can 
be traced to the contribution of the field associated to the phase of the last component 
of the field z^, which cannot be eliminated by a gauge transformation when x / 0. In 
turn the difference between Eq. ( |14.16|) and Eq. (|14.17|) is originated by the unsuppressed 
contribution of the tadpole graphs involving closed loops of vector propagators, shown in 
Fig. Trivial power counting arguments show that no inverse powers of 1 + x/ appear in 
the perturbative evaluation of these diagrams. 

The agreement between perturbative and 1/N evaluation of the ratios of A parameters 
is a strong confirmation of the commutativity of the 1/ f3 and 1/N expansion with respect to 
renormalization-group properties of the models. Even the apparent singularity of the x — > 
and x — > 00 limits has no consequences on the exchange of the (3 — > 00 and iV — > 00 limits. 
Further confirmations can be obtained by the comparison of the perturbative /3-function 
coefficients with those obtained by expanding the resummed 1/N lattice (3 function, which 
can be easily obtained from 8m\. 
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B. Lattice (3 and 7 functions 



Let us come to the evaluation of the 0(1/ N) lattice renormalization-group (3 function. 
We may apply the homogeneous renormalization-group equations to the expression of the 
mass gap 



m 2 (M,f) = 0. 



(14.22) 



We can expand 13(f) in powers of 1/N in the form 

Our choice of lattice action allows us to use the relationship 



(14.23) 



M l exp ( -j- ] . 



(14.24) 



in Eq. (|14.22|) to obtain 
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Further substitutions in Eq. ( |14.22j ) lead to 



/w> = w/)]1!g 



(14.25) 



14.26) 



and, since we know f3\ in the SM scheme, we immediately obtain 



?(SM) 



P p d (8m{ 



7r 2n df \ m\ 



(14.27) 



Eq. (|14.27|) admits a natural interpretation. We must recognize that a change in the regular- 
ization scheme corresponds to a reparametrization of the model, i.e. /' = /'(/)• Covariance 
of the renormalization-group equations under reparametrization implies 



W) = P(f(f)) 



df_ 
df> 



-1 



14.28) 



As a consequence Eq. ( |14.27|) implies 

1 



./ 



(SM) 



N 2 Wq \J\I z 



(14.29) 



The improper integral obtained by substituting Eq. ( |1 1 . 52| ) into Eq. ( p.4.27 ) is defined 
according to the prescription ( |11.70|) . It is worth noticing that all the residues in the complex 
integration vanish, in contrast with Eq. ( |11.52| ) itself. Taking the derivative with respect to 
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/ = 1/(2/3) in Eq. ( |1 1 . 52|) is completely straightforward and we shall not write down the 
result in the most general case. We shall however consider a few interesting special cases. 
For x = (CP^ -1 models) we have 
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For x^O and C2 = we have 
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For x — > oo we have 



7T 



AT 



AT 2. iVi (2^2 ^a) + ^w 



1 

+ TT 



d 2 fc 1 x ^ 
9 2^ 



1_ f M 'dk 2 lrr 



N J (2tt) 2 2 ^ (CV + /3A) M ) 2 AT,/ 4vr fc 2 ln 2 (£; 2 /m 2 



(14.32) 



Finally, for 0(2A^) models we have 
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P(f) 
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The evaluation of the lattice renormalization-group function 7 follows essentially the 
same path. From Eq. ( |5.41 ) we obtain the relationships 
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Eq. ( 14.35 ) in turn is consistent with 







where 



7 W (/) =7 (SM) (/ (SM) (/)) _ p{f) ^_ lnm 
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(14.34a) 
(14.34b) 



(14.35) 



(14.36) 



14.37) 



is the additional finite field-amplitude renormalization due to the change of regularization 
scheme. We may appreciate the fact that Ci(/) = 0(f 2 ), as expected. 

15. FINITE SIZE SCALING 



The study of finite-size effects is quite important, both from a purely theoretical point of 
view and in the context of controlling systematic deviations from the infinite-volume limit 
in numerical simulations. Finite size scaling in the large-iV limit has been widely studied 
for different geometries and space dimensionalities by Brezin and collaborators ||155| , |156 



A finite-volume approach was applied to the problem of evaluating the low-lying spectrum 
in two-dimensional spin models by Luscher ||157|| and extended by Floratos and coworkers 
158- 1 . The systematic analysis of 1/JV finite-size effects is however rather recent [161 



Let us review the main results of this analysis. 
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Any coordinate-independent physical quantity Q defined in the context of the 1/N- 
expandable finite-lattice model will in general depend on four different parameters: 



Q = Q(f,a,L,N), (15.1) 

where L d is the physical volume in d dimensions and a ~ 1/Mt, is the lattice spacing. In 
the infinite-volume limit and in the scaling region (according to the discussion presented 
in Sect. |14]) all separate dependence on / and a can be made disappear by parametrizing 
everything in terms of the physical mass gap m 2 (a, f, N), solution of Eq. ( |14.22| ). The finite- 
size-scaling relation stems from the observation that one can reach the infinite-volume limit 
(L/a — > oo) while simultaneously keeping a constant finite value of mL. As a consequence 

Q(f,a,L,N) {Q) 

— — > f w >(mL,N) 15.2 

Q{f, a > oo, N) /-*o 

mL=const 



The 1/N expandability in turn implies that, assuming 



Q(f, a, L, N) = Q (f, a,L) + ^ Q 1 (f, a , L) + O (-^ ) (15.3) 



and 

m(f N\ = mJf\ 4- _ ,,, w t ^ 

we may expand the finite-size functions f^> in the form 



m(f, N) = mo{f) + ^- T m 1 (f) + 0[-^), (15.4) 



fW(mL,N) = f(mL) + i f[ Q \mL) + O (J^J . (15.5) 
Substituting Eqs. QT5T3; ), ( |T5^D , and (|T5T5| ) into Eq. flloT^ ), we obtain 

Q {f, a, oo) 

, a, oo) _ ( mQ L) ( r I}± \ /o (Q) ( m o-^) 



fi Q) {m L) Q {f,a,L) Q Q {f, a, oo) 1 \m J f^ Q \ mQ L) 

Eq. ( |15.6| ) is the most general form of the 1/iV-expanded finite-size-scaling relation. In 
order to gain further insight, one must consider the specific properties of the quantity under 
investigation. In any case, a basic tool in the analysis is the knowledge of the finite-size- 
scaling properties of the finite-size mass parameter defined by the gap equation: 

Ly 1 _ 1 f d 2 g 1 , , 

£2 2^2 + m 2 2f J ^2 q2 + m 2^ ■ 

where the sum runs over the momentum lattice modes, i.e. = 0, 2tc/L, 2n(L— 1)/L. 
In the case c 2 = and in the scaling region, defining z L = m L L and z = m L, we could 
establish the relationship 
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z = z c exp {~u)(z L )} , 



(15.8) 



where z c = 4.163948 and in the region zl < 27r the function to may be defined by 
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4 OO 
4-7T 

(z L ) = — + A7cJ2(- l ) n 4 n d n+ i, (15.9a) 



n n=l 



dn = ~, rr; — ^5 kt~ , n > 1. (15.9b) 

y ' ni,n 2 =-oo v 1 Al 
(ni,n 2 )^(0,0) 

The function Zq(zl) is monotonic and invertible. Therefore all subsequent calculations can 
be performed making use of the auxiliary variable zl, which simplifies many computations. 

Without giving further technical details, we mention that finite-size-scaling functions 
were computed in the 1/N expansion of masses and magnetic susceptibilities, both in the 
pure O(iV) case [161] and in the pure CP 7 ^ 1 models (JT62J. In O(N) models the 1/N ex- 



pansion of finite-size functions is an accurate description of finite-size effects in all possible 
regimes: small volume (mi <C 2tt), where the results can also be compared, by asymp- 
totic freedom, with those obtained from finite-volume weak-coupling perturbation theory 
]; large volume (mi 3> 2n) where, due to the existence of a physical mass gap, one 



expects exponentially fast convergence to the infinite- volume limit ||166|| ; and intermediate 
volume. 

In CP^ 1 models at intermediate volume and for N not too large (N < 100), new 
phenomena occur: not every physical quantity is expandable in a 1/N series; moreover, 
even if we limit ourselves to 1/iV-expandable objects, we must observe that the scale of 
finite-size effects is not set by the mass gap, but instead, since we are in presence of a 
"weak" confining potential, it depends on the (semiclassical) radius of the bound states; the 
radius in turn grows like iV 1//3 and therefore does not have an analytic dependence on 1/N, 
a result not unexpected in the light of the semiclassical results on the bound state spectrum 
presented in Appendix [E[ 

Another subtle point in the study of finite-size effects in CP^ -1 models is related to the 
properties of the Abelian Wilson loop on finite lattices. We only mention here that, defining 
the Polyakov ratio, corresponding to the derivative of the static potential introduced in 
Sect. |, 

in the infinite-volume limit and in the scaling region one should find the Abelian string 
tension 

xp(R) — %— -0-^. (15.11) 
However, on finite lattices, even in the scaling region, 

lim Xp(R,L,a,f,N) = /(P) ^ ^ {m ^ R/L) (15 u) 

mo^O O N—>co 
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and one may show that 



lim fP(m L,R/L) = l-^. (15.13) 
We may define a function measuring the deviations from the infinite (periodic) volume limit: 

g < f , = gwi . (1514) 

(?q in turn can be shown numerically to enjoy a factorization property: for large m^L and 
mo-R 

(m L, R/L) « 1 + 0(woL) ip(R/L), (15.15) 

which can be understood in terms of an effective Yukawa interaction replacing the Coulomb 
potential on finite lattices. 

If we compare Eq. ( 5.58 ) with the definition ( [13.14 ) of the topological susceptibility, we 



recognize that the latter quantity is strictly related in CP^ 1 models to the Abelian string 
tension, i.e. 

(cf. Eq. (|5.72|) ). Therefore, as a side effect of the above analysis, we are led to the (not unex- 
pected) result that Xt should vanish on any finite lattice. However, since the infinite-volume 
limit is reached smoothly, it should always be possible to devise an appropriate limiting 
procedure to extract infinite- volume information from finite- volume, finite-a measurements. 
Effects of topology in finite volumes were also studied, for different geometries, in Ref. ||167|| . 

In closing the present Section, it is relevant to observe that finite size scaling in 0(A) 
models has been the subject of studies concerning the three-dimensional case, where a 
second-order phase transition occurs at a finite value (3 C of the coupling. In Ref. [|168|1 finite 



three-dimensional lattices were studied in the context of the 1/N expansion. The method for 
treating near-critical behaviors in three dimensions, originally developed in Ref. ||169fl for the 
infinite- volume limit, is very reminiscent of the asymptotic expansion techniques employed 
in the present work. 

Another approach to finite-size effects in three-dimensional 0(A) models near criticality 
was developed by Hasenfratz and Leutwyler employing the techniques of chiral perturbation 
theory 



16. HIGHER ORDERS OF THE 1/N EXPANSION ON THE LATTICE 

We have not seriously addressed the problem of evaluating 0(1/N 2 ) contributions in 
the scaling region. A finite-/?, finite-lattice calculation can certainly be performed with no 
conceptual problems, but with some technical troubles in evaluating accurately two-loop 
lattice integrals involving dressed propagators. 
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This approach has been put forward in recent years by Flyvbjerg and collaborators 
171H173|| , who explicitly studied the case of O(N) models. They evaluate the two-point 



function of the O(iV) nonlinear a models up to 0(1/N 2 ), on finite square lattices and 
for fixed values of N, typically N = 3,4, in order to compare with existing Monte Carlo 
simulations. From the two-point function they can extract the numerical value of such 
physical quantities as the mass gap and the magnetic susceptibility. The comparison with 
high-precision Monte Carlo results allows an estimate of the systematic errors involved in the 
series truncation to zeroth, first, and second order. These errors appear to be uniform and 
smaller than the expected magnitude of the neglected terms. Further insight is obtained by 
the use of Fourier-accelerated numerical evaluation of Feynman diagrams and extrapolations 
of finite- volume results to infinite volume by phenomenological finite size scaling |174|J175| . 

These improved results lead to agreement with Monte Carlo data, within the expected 
errors, for N > 3, and give for N = 3,4 extrapolated mass gap - yl-parameter ratios 
consistent with the exact continuum results of Ref. | 112 |, reported in Eq. ( 4.15| ). 

However a fully analytic approach to higher orders in the 1/N expansion would require 
extracting the scaling contributions along the lines defined in principle in Sect. [TTJ. This 
extraction in turn would involve a proper treatment of the regularization problem, which 
may not be straightforward in presence of higher loops, if we want to stick to our favorite SM 
scheme: we might run into technical problems similar to those involved in generalizing the 
BPHZ scheme to massless theories. In any case, based on the proofs of renormalizability of 
the 1/N expansion, we believe there should be no general obstruction to such a calculation. 



17. A DIFFERENT APPROACH: SCHWINGER DYSON EQUATIONS 

The approach to large N based on the effective action and effective Feynman rules is 
by no means the only way of generating an expansion that is naturally organized in powers 
of 1/N. Writing down Schwinger-Dyson equations for U(iV) (0(2iV)) invariant correlation 
functions, it is possible to recognize that N occurs only polynomially in the coefficients of 
the equations themselves. It is therefore possible to truncate the (a priori infinite) set of 
Schwinger-Dyson equations by keeping only terms and equations down to a chosen power 
of N. When the resulting finite set of equations is solved, the solution depends on 1/N 
through all powers, and it is equal to the sum of an infinite subseries of the exact 1/N series. 
It is therefore at least as accurate as the corresponding truncated 1/N series; in practice one 
can get sensibly higher accuracy, as shown in the original papers by Drouffe and Flyvbjerg, 
explicitly concerned with 0(N) models ||165| , |176| , |177|| . 



The derivation of the Schwinger-Dyson equations is essentially straightforward in the 



generating functional formalism, and we refer to [177] for all details. We just present in 



Fig. |24| the graphical form the equations involving the fundamental field and the a-field 
propagators in the 0(2iV) model, truncated to 0(1/ N 2 ). There is an ambiguity in the 
location of the bare and dressed vertices; we fixed it by the prescription that each a-field 
propagator should connect a bare and a dressed vertex. 

We tested our prescription by a very simple and illuminating example, the exactly- 
solvable continuum one-dimensional 0(2iV) model. Since the 0(1/ N) truncation amounts 
to identifying the bare and dressed vertex, we may try to solve the resulting equations by 
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the Ansatz 



G-\p) = A{p 2 + m 2 ). (17.1) 

Simple integrations (performed for convenience in dimensional regularization) lead to an 
explicit solution in the form 

m = W' (17 - 2a) 

2/V 

G-\p) = w ^( P 2 + m 2 ), (17.2b) 
2N — 1\ 2 N 



^ 1(P) =(^J Tn^T^- (17 ' 2C) 

Actually Eqs. ( |17.2| a) and (|17.2[ b) correspond to the exact solution of the one-dimensional 
0(2iV) model ||178|| , thus showing the power of the approach: albeit truncated to 0(1/N), 
the equations enable us to resum the whole 1/N series. 

We must now verify that the 0(1/N 2 ) truncation does not spoil the result. However, 
because of our prescription, this is simply achieved by the Ansatz 

N f2N-V 2 



A- 1 (p) = — [ — — -y(p), (17.3a) 

KFJ m\2NJ p 2 + 4m 2 yFh v ; 

V ia) (p,p') = V(p-p'), (17.3b) 
insuring that the value of G~ l (p) is unchanged, while by consistency one can find 



2N V P 2 + Im 2 

Eq. ( |17.4|) can be verified in the context of the standard 1/N expansion. 

Although rather promising, the Schwinger-Dyson approach has till now only be applied 
to the 0(1/N) truncation of the two-dimensional O(N) models ||177| ; by comparing their 
results with high-precision Monte Carlo data for the 0(3) and 0(4) models, the authors find 
an apparent uniform systematic error of 0(1 /N 3 ). Higher orders and more general mod- 
els involve, besides the abovementioned problem of locating the dressed vertices, technical 
difficulties related to performing the higher-loop integration with the needed accuracy. 



18. FERMIONIC MODELS 



The next obvious extension of the 1/N approach in the study of model field theories is 
considering fermionic degrees of freedom. The natural counterparts of the bosonic models 
we have discussed are theories with four-Fermi interactions and U(N) symmetry, power- 
counting renormalizable in two dimensions and 1/N expandable. These theories are ultra- 
violet renormalizable order by order in the 1/N expansion in less than four dimensions, as 
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shown by Rosenstein and coworkers ||1 79 - 181 1 and by Kikukawa and Yamawaki ||182|| , essen- 
tially because their ultraviolet limit has the same relevant operator content as the infrared 
limit of the superrenormalizable Yukawa model. These arguments were further developed in 
Refs. Ul83| - |185|l , and the whole subject is carefully reviewed in the introduction of Ref. ||186| . 
Critical indices were computed by Gracey to 0(1/N 2 ) ||187| . 

A sufficiently general U(iV)-invariant two-dimensional continuum Euclidean Lagrangian 
depends on three couplings | |188| |: 



£ = ffitf, - lg s (^) 2 - lg p (^75^) 2 - \g v (ip^^Y 



;is.i) 



where ip is an iV-plet of Dirac fermions. This Lagrangian interpolates between the 0(2iV)- 
symmetric Gross-Neveu model (g p = g v = 0) and the SU(iV)-symmetric Gross-Neveu model 
(9s = 9p = Ng v ), enjoying a global axial U(l) invariance. 

Gross-Neveu and chiral Gross-Neveu models possess factorized S'-matrices in two dimen- 
sions, and other exact results can be obtained, in full analogy with their bosonic partners. 
The related literature has been steadily growing in the last twenty years, and we shall not 
even try to give references for the continuum formulation, addressing the interested reader to 
the (partial) bibliography appearing in Ref. [189]. However, as already observed for bosonic 



models, only a minor effort was done in going beyond the leading large- N approximation, 
both in the continuum and in the lattice versions of the models. 

Lattice formulations, as it is well known, are plagued by the supplementary problem of 
fermion doubling. The problem is solved in principle in any dimension by the introduction 
of the Wilson term [ |190|| , which however leads to an unavoidable complication in both 
analytical and numerical computations. In the present context we only want to summarize 
the available results concerning the large- N limit and the 1/N expansion of lattice Gross- 
Neveu and the chiral Gross-Neveu models, and include a new result coming as a rather 
natural extension of our previous analysis. 

The first lattice formulations of the Gross-Neveu models admitting the correct contin- 

|19T 



192 



therefore 
fermionic 



and involved staggered fermions 

naive' 

A lattice formulation with Wilson fermions of the chiral 
Gross-Neveu model was introduced and discussed in the large- N limit in Ref. [|194j[ (see also 
Ref. ||195|| ). The Gross-Neveu model in the Wilson formulation was studied to 0(1 /N) in 
the seminal paper by David and Hamber [196 , where the notion of an asymptotic expansion 



uum limit were presented in Ref. 
they described models with 0(4iV) symmetry, N being the number of 
components (cf. also Ref. [|193|| ). 



of the effective propagator was introduced. In Refs. [|197| , |198|| the Symanzik improvement 
program was applied to the large- N Gross-Neveu model with Wilson fermions. A systematic 
analysis of the 1/N contributions to the lattice Gross-Neveu model was finally performed 
in the staggered version in Ref. ||199|| , and in the Wilson-Symanzik version in Ref. ||200| . 
For completeness we must mention that a lattice formulation of CP^ -1 models coupled to 



fermions was discussed in Ref. [[201 ], and that the problem of formulating lattice versions of 
two-dimensional 1/iV-expandable supersymmetric models was addressed by a few authors 
134|,|202||, but never systematically investigated. 



Our new result concerns the possibility of obtaining an integral representation of the 



effective propagator A 
action 



in a staggered version of the Gross-Neveu model, defined by the 
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x,y 



2 



;is.2) 



where \ x is a iV-plet of (one-component) fermionic fields, D xy is the Susskind "differential" 
operator 



and the fermions are coupled to the Lagrange-multiplier field o by 

In this case, we can obtain an integral representation of the fermionic integral 



;i8.3) 



;is.4) 



A7\ 
(°0 



tr 



I 



f 



d 2 p 

(2tt) 2 I j^p+i^ _|_ mo i^^p-^k^ + m 
d 2 p 



(18.5) 



where p M = sinp M , and the mass parameter mo is momentum-independent. By repeating the 
arguments of Sect. [|, we find that ATs can be computed in closed form along the principal 
diagonal of the momentum lattice: 



ir 1 + m 2 



mf, — cos I 1 
+ 



\ ""0 

4 1 1 



ml + l 



n 



cos 2 / 



cos/; V(l + m2) 2 ' l + m 2 



K 



7T 1 + TOq COS Z \ 1 + m Q 



;is.6) 



More generally, using the standard Feynman parameter representation we obtain 

d 2 g c + V b^cosq^ 



1 + m 2 , - <fy cos g M 



2 ' 



(18.7) 



where 



2 

COS fc 



4a, 



c = m r 



By essentially trivial algebraic manipulation we are led to the final result 



(18.8a) 
18.8b) 
(18.8c) 
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_ ml + c(l + 2ml) f 1 , 1 C 3 ^(C) 
<"> W " 4vr 7o («i« 2 ) 3/2 1 - C 2 



i /- 1 . c 6„ r s(o 



2tt J Q aX (a 1 a 2 )V2 2-, « 

where 



A 4 P 



i-C 5 



+ E(() - 2K{Q 



;i8.9) 



4aia2 



(l +m 2)2_(a 1 _Q 2 )2 



;i8.io) 



19. CONCLUSIONS AND OUTLOOK 

In our opinion the most important conclusions that can be drawn from our results are 
summarized by the following statements. 

1. Nontrivial asymptotically- free two-dimensional Euclidean field theories can be con- 
structed, in the context of the 1/A expansion, starting from a lattice formulation and 
exhibiting explicitly the existence of a scaling region. The accuracy of our construction 
is 0(1/N), but there is no obstruction to higher-order extensions. 

2. In the scaling region, results that are expressible in terms of adimensional ratios of 
physical quantities are universal, i.e. they do not depend on the specific lattice model 
adopted as long as the physical parameters are kept fixed. Moreover these results are 
unaffected by the pathologies of standard perturbation theory and can be unambigu- 
ously predicted. 

3. The width of the scaling region however necessarily depends on the choice of a lattice 
action. In turn, it is widely independent of the 1/A corrections, since these depend on 
the effective propagators and vertices, whose scaling properties are fixed by the (large- 
N) effective action and are modeled upon the scaling properties of the large- N lattice 
mass gap (cf. Eq. (|11.8 ) and Fig. [15]). For standard nearest-neighbor interactions, 



scaling within 10~ 3 is achieved starting from 2f ~ 1.25 (J3 ~ 0.8), corresponding (for 
not too small N) to a correlation length l/(ma) ~ 27. 

4. "There ain't nothing like asymptotic scaling" in the real world (excluding very large 
N). The asymptotic scaling region is in our language the small-/ region, where the 
behavior of mf/ml is well approximated by the two-loop perturbative renormalization 
group, i.e. its finite part (in the notation of Sect. ||) is very close to its value at / = 0. 
As an example, for the 0(A) models (the "best case") at c 2 = 0, with A as large as 
20, the mass gap is approximated by the (two-loop) asymptotic formula within 1CT 3 
for /3 > 3.3, i.e. I /{ma) > 10 s . 

5. Perturbation theory may however be a good guide to the physics of the models, in 
that it commutes order by order with the 1/A expansion and it leads to the same 
renormalization-group functions and asymptotic behaviors. Moreover, by summing 
over a sufficient number of perturbative terms one may reproduce the correct lattice 
A parameter, renormalization constants and perturbative tails throughout the whole 



scaling region (cf. Figs. p!8|-[2T|) 
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As we have shown, the 1/N approach may be successfully extended in many different 
directions. The major limitation we could not however bypass is the restriction to models 
where the fields belong to the fundamental (vector) representation of the symmetry group. 
The problem of extension to fields in the adjoint (matrix) representation, like principal chiral 
models and gauge theories, has been the stumbling block of the 1/N expansion in the last 
decade. A breakthrough in this domain could turn the 1/N expansion from a toy in the 
theoretical playground into a major tool in the analysis of realistic physical models of the 
fundamental interactions. 
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APPENDIX A: PERTURB ATIVE RESULTS 



In the literature on perturbative calculations, it is usual to report the results in terms 
of a rescaled renormalized coupling t, following the notation first adopted in Ref. [II]. All 
four-loop-order (3 functions of nonlinear a models on symmetric spaces can be found in 
Ref. |0J. Here we are only interested in two special cases: 

1. 0(A)/0(A-1) spaces 

P{t) =et-(N- 2)t 2 [l+t + \(N + 2)t 2 + (-^A 2 + fN - f + |(A - 3)((3)) t 3 ] . 

(Al) 

The result for 0(A) models is obtained by setting t = l/(27rA/3 v ). 

2. U(A)/(U(A-l)xU(l)) spaces 

f3(t) 2*et- Nt 2 [l + 2t+ (|A + 2) t 2 + (§A 2 + f N + l) t 3 ] . (A2) 

The result for CP^ -1 models is obtained by setting t = l/(2irN(3 g ). 
The local nonderivative scaling operators can be expressed in terms of orthogonal poly- 
nomials in the variable a 2 = Z\Z\ ||0| (see also |[122|| ). Their anomalous dimensions were 
computed to four-loop order: 

1. 0(A) models 

(it -1) 

The scaling operators are the Gegenbauer polynomials C h 2 (a) and 

7l (t) ^/(iV + /-2){t+|(A-2)t 3 

+ (A - 2) [-IA + | + |C(3) (1 - ll(N + I - 2))] t 4 }- (A3) 

2. CP^ 1 models 

The scaling operators are the Jacobi polynomials Pjf* 2,0 '(2ct 2 — 1) and 

7 fc (t) = 2A;(A + k - 1) {t + |At 3 + (A + 6) [§A + |((3) (A — k(N + k — 1))] t 4 } . 

(A4) 
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The e-expansion of the critical exponents may be extracted from the above results by 
finding the critical point t* defined by f3(t*) — and applying the relationships 



77 = -e + 7i(0, 
1 



v 



P'(t*) ' 



(A5a) 
(A5b) 



holding for O(N) models. 



APPENDIX B: EFFECTIVE PROPAGATORS IN d DIMENSIONS 

The mass-gap equation takes the form 

_ r (l-\d) 2^/2-1 

The inverse propagator of the a field is 

d d p 1 1 



(Bl) 



(2n) d q 2 +ml (p + g) 2 + m 2 , 

r(2-|d) r 1 dx 



(4vr) d / 2 A \p 2 x(l- 



( 47r )rf/2 



[p 2 a;(l-:r)+m 2 ] 2 ^ /2 

1 , 

2 '2' 2' e 



r ( 2 -H ,1-2 . 1,13 I 



(B2) 



where F is the hypergeometric function. For p — 0, the inverse propagator assumes the 
value 



^-1(0)- r ( 2 ~H f m 2^/2-2 

Special values for the lowest integer dimensions are 



(B3) 



^ = 0: A^{p) 



d = 


1 : 




d = 


2 : 




d = 


3 : 




d = 


4 : 





8 1 

We 
i 



1 + lm^ + 1 



2/5 



2/5 



+ 



4£ 



In 



p 2 + 4m 2 ( p 2 + 4m 2)2 

(B4a) 



mo p 2 + 4m 2 , p 2 + 4m 2 , ' 



1 1 h f +1 



1 



C ln C + l 

2np 2 £ £ - 1 p 2 + 4m 2 2vr 11 ^ - 1 ' 

1 p 
arctan 



47rp 



2mo ' 



(4tt) 2 



C - 1 



(B4b) 
(B4c) 
(B4d) 
(B4e) 
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The inverse propagator of the 6 field is 



where 



d d p q ■ (p + q) — m. 



= 2 



= 2 



m; 



2\d/2-l 



d d p 

(27r) d q 2 + m 2 d-lj (2n) d [q 2 + m 2 ) [(p + g) 2 + m 2 
(47r) d / 2 

r(i-irf) 

{Any/ 2 



o [p 2 :r(l — x) + ml] 1 d ^ 2 



1,13 1 



For p = 0, we have 



(B5) 



(B6) 



^(0) = 0. 



(B7) 



APPENDIX C: CONTINUUM INTEGRALS 

A number of continuum integrals occurring in the evaluation of the constants appearing 
in 0(1/N) results can be computed analytically. A typical dimensionless SM-regulated one- 
loop integral in the 1/N expansion, after angular integration has been performed, takes the 
form 



J-s = m ~ a- 2 



(CI) 



where T^ UV ^F are the terms in the series expansion of F in powers of m\ that are only 
ultraviolet-divergent ( "perturbative tails"), while T°F are the terms that are both ultraviolet 
and infrared divergent; A is the canonical dimension of F. F{p 2 ) m'^) may be represented in 
the form 

m o P (^)= m om, (C2) 

where £ = a/1 + Am\jp 2 . 

It is now convenient to perform the following changes of integration variable: in the first 
integral we set 

t= J=V dp 2 = mg (l - df; (C3) 
in the last two integrals we set 
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u= f_ dp 2 = m 2 d (C4) 

ml t F v 1 



By observing that t — u + 2 + 0(l/u) , we obtain the representation 



^4+2 / A / f+1 N M 



47r/ reg = lim ^ / dt 1 - - / 







dM T ( uv )F(m) - [ A du T°F(u) \ . 



(C5) 



Now, in order to parametrize explicitly the regulator-dependent part of J reg we can split 
the last integral into the two regions M 2 jm\ < u < exp(a) and exp(a) < u, with a = 1 for 
integrals involving and a = 3 — 27rx for integrals involving The integration 



f-cxp(a) 

is now trivial, and we are left with the task of evaluating 



/•exp(a) 

/ duT°F(u) (C6) 

J M 2 /mi 



(C7) 



Exact analytic results have been obtained in the following instances: 



/(0=lnln 



71+2 / 1 \ l" 71 fdu 2 

c = lim <! / dt I 1 ) In lnt — / dwln \\nu' 



/(0 = 



yl-»oo y t 2 / Je U l nM 

/■°° dt Z" 1 Z" 00 dt 

— / —lnlnt — / dw In |lnu| = — 2 / — In lnt 

2ie; (C8) 



111 



lim i rv^ * 2 1 r * 



yi->oo y t 2 / lnt Jg m In if J y 1 t 2 

2je; (C9) 
1 4£ 
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; A + 2 ( t - if 3 1 ^ dw 3 

c = lim ^ - / dt ± — — ■ + / — 



yi^oo [ J 1 t hit t 2 — t + 1 J e u \nu 

r°° 1 - t 2 

-3/ dt ? mint = -3( lE - Cl ), (CIO) 

A (t 2 -t + i) 2 u J K J 



where c\ is given by Eq. (|5.38|) ; 



/(0 = - v e 



21n.^ 1 ' 



^-1 

f Z^ 2 , t - 1 1 f^du 1 
c = iim < / dt -=r -f 

Ui tv^ lnt 7o 



lim < — 



_ m ln« 
00 dt 1 /' 1 " £ dw 1 



0. (Cll) 



APPENDIX D: EFFECTIVE VERTICES IN THE CONTINUUM 



The effective vertices of the 1/N expansion are nothing but one- loop integrals over the 
fundamental field propagators with appropriate couplings to the external lines. The problem 
of evaluating the most general continuum one-loop integral in two dimension is solved in 
principle in terms of elementary functions [|203| - |205l . It is however convenient to derive ex- 
plicit expressions for those special kinematic configurations entering the actual computations 
we would like to perform . 



One basic ingredient is the three-point scalar vertex 



V 3 (pi,P2) 



1 



d 2 q 1 

(2n) 2 q 2 + ml (q + pi) 2 + ml (q + p 2 ) 2 + m 2 . 



2 ■ 



(Dl) 



a symmetric function of pi, p 2 and p% — p 2 . Two-dimensional identities allow a reduction of 
the integrand to a combination of terms involving only two fundamental field propagators. 
The integration is then straightforward, and the result is 

+ pl{Pi{Pi ~ P2))A-^(p 2 ) + (pi - p 2 f{p 1 p 2 )A-^(p 1 - p 2 )] , (D2) 

where 

D(pup 2 ) = pIpKpi - P2) 2 + 4mjj \p\p\ - (P1P2) 2 ] • (D3) 

Let us now consider the four-point vertices: the exceptional configurations we are inter- 
ested in are the cases when the external momenta are equal two by two. Let us define 
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d 2 g 1 1 1 

(2vr) 2 [g 2 + mg] 2 (q + Pi) 2 + m 2 (g + p 2 ) 2 + m 



(D4) 



Again by applying algebraic identities we are led to an explicitly integrable expression. The 
final result is 



v} a \ Pl ,p 2 ) = D- l ( Pl ,p 2 ) 



(pj - (pip 2 ))pj 
p 2 + 4ml 

„2 f m m \ \ „2 



<)(Pi) + ^(0) 



+ 



+ D- 2 (p l ,p 2 ){[{pi -P2) 2 (pm) +p\pI - (pm) 2 ] 

x (pt - (P1P2)) p? ^ (pi) + (p? - (P1P2)) P 2 

+ [(Pi -P2) 2 (P1P2)] 2 ^J(pi -P2)] 

- D- 2 ( Pl ,p 2 ) [ V \v\ - ( Pl p 2 ) 2 ] Up\ + 4m 2 1 )((p 1 p 2 ) - p?)^(pi) 



+ (p 2 , + 4mg)(( Pl p 2 ) - p^)^ a 1 ) (p 2 ) 

+ [(pi - P2) 2 + 4mg] (pi - p 2 ) 2 ^(^(pi - p 2 )}- 



We must also evaluate 
d 2 g 



^f(Pi,P2) 



1 



(27r) 2 q 2 + ml (q + pi) 2 + ml (q + p 2 ) 2 + ml (g + pi + p 2 ) 2 + m 2 



One can show that the result is expressible in the form 



vi b \ Pl , P2 ) 



(P1P2) 



[V 3 (pi,p 2 ) - V 3 (pi, -p 2 )] 



(D5) 



(D6) 



(D7) 



In order to compute the correlation function of the composite operator Pij(x), we also 
need the mixed four point scalar-vector vertices in exceptional momentum configurations. 
We quote here the definitions: 



d 2 g 



{2q IM + k )1 )(2q u + k u ) 



(2tt) 2 (g 2 + m 2 ) 2 [(g + p) 2 + m 2 ] [(g + k) 2 + ml] ' 
d 2 g 1 2q tl + k tl 1 2g^ + 2p u + k v 



(27r) 2 g 2 + m 2 , (g + k) 2 + ml (g + p) 2 + ml (g + p + k) 2 + m\ 



(D8a) 
(D8b) 



Actually we only need the combination of vertices appearing in Fig. [7] and this can be 
shown to be a transverse tensor. Therefore we can limit ourselves to computing 



kfi k v 



) V^(P, k) = -{k 2 + 4m 2 )\/ 4 (a) (p, k) + 2V 3 (p, k)+(l- 

(p-k) 



Am 2 (p-k)\ 



A- 2 



P 



p 2 + 4m 2 , 



(p + /c) 2 + 4m 2 1 + 



A7\(p) 1 



F(p 2 + 4m 2 ) k 2 W 



(D9a) 
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and 

(V - kj ^) Kt\p, k) = -{k 2 + 2p 2 + 4n$)V<%, k) - 1 A^(p) + 1 ^(p - fc) 

+ 1 Z\^(p + k) + 2^ 3 (P, fe) + 2V 3 (p, -k) . (D9b) 



APPENDIX E: THE BOUND STATE EQUATION IN THE LARGE- AT LIMIT 

As discussed in Sect. |5], we must solve the following eigenvalue Schrodinger equation: 

-^ + ^Mx)(l-e~ xp )^ = ^ (El) 

where p = m R and x = m g /m . The resulting bound state masses will be = m (2 + e). 

It is convenient to introduce new "natural" variables: a rescaled coordinate y = xp, a 
"weak" coupling 

r \ V 3 /a \ V 3 



»=ljv^W ^ hvl • (E2) 



and a rescaled energy eigenvalue rj = e/g 2 . Eq. ( |El| ) now turns into 

d 2 ^ 



dy 2 



+ t 3 (l-e~ y )ip = 7]t 2 ilj, (E3) 



where t = g/x = grrio/mg. 

The general (unnormalized) solution of Eq. (|E3|) is a Bessel function: 

Mt,y) = J2 tV i^{2t 3/2 e-y/ 2 ), (E4) 

and the eigenvalue condition simply amounts to 

J 2tvt^ ( 2t3/2 ) = (even-parity levels), (E5a) 
^v^(2t 3/2 )=0 (odd-parity levels). (E5b) 

More specifically, denoting by the fcth zero of J„(z) and by j' vk the fcth zero of J' u (z), 
the nth energy level r] n is determined by solving 

jW^,(n + i)/ 2 = 2t 3/2 (oddn), (E6a) 
J 2tv ^r, n/2 = 2t 3/2 (even n) . (E6b) 

The rescaled energies r\ n are plotted as functions of t in Fig. |25|. The energy eigenvalues e n 
themselves and their dependence on the mass parameters can be easily derived from r) n (t). 

This kind of analysis is especially suitable for discussing the large-t behavior, correspond- 
ing to the condition 
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6tt^ 1 ' 



m B < gm = ( — J m . (E7) 



From the asymptotic formula for j[ 1 at large order Z [pOC 



one obtains an expansion in the form 



k=i 



(E8) 



fc=0 

where 70 = 1.01879297. Similar expansions can be derived for higher energy levels. As long 
as t 3> 1 these are good descriptions of the mass spectrum, and in particular when t — > 00 
they reproduce known results for the bound states of the CP^ -1 models. In passing we 
notice that the condition ( |E7p can be rephrased in the form 

m e (R) B < 1, (E10) 

where (R) B is the semiclassical bound state radius. Eq. ( EIQ ) is an obvious consistency 
condition for the calculations. 

Eqs. ( |E5| ) indicate that, for sufficiently small values of t, higher bound states may disap- 
pear from the spectrum. More precisely, defining t n by the condition 

Jo,(n + i)/ 2 = 24 /2 (oddn), (Ella) 

Jo,n/ 2 = 24 /2 (evenn), (El lb) 

the corresponding state satisfying r\ n = t n , the nth level disappears for all t < t n . In 
particular, when t < t 2 — 1.130756402 all excited bound states disappear. This phenomenon 
is illustrated in Fig. |25|. 

We may appreciate that, as long as m# ~ grriQ, one may show that the mass of the 
fundamental particle is well approximated by 

m% = ml + ^ E x {-ml) = ml + I ^ , (E 12) 

implying 

— = l + -9 2 t. E13 
m 2 

The threshold condition for the nth bound state in turn has the form 

^ = 2 + g\ n = 2 + gH n , (E14) 
m 

and therefore it amounts to the condition 

m B (n)(t n ) = 2m F (t n ), (E15) 

as one may expect on physical grounds. 

Finally, the threshold condition for the first bound state is simply t = 0, which corre- 
sponds to g = 0, and from Eqs. ( |E2|) and ( |5.63| ) this implies me = 2m , a result known by 
independent arguments. 
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APPENDIX F: LATTICE INTEGRALS 



We know from general theorems that the perturbative expectation values of quantities 
invariant under the full symmetry groups of the models must be infrared-finite. Therefore 
in principle it must be possible to represent finite lattice expectation values in terms of 
lattice integrals only. However our regularization technique has led to the introduction 
of continuum counterterms. There must therefore exist (infinitely many) simple identities 
connecting lattice integrals and their continuum counterparts. 

We have found all the identities that might be relevant in a three-loop perturbative 
computation of lattice renormalizat ion-group functions (when C2 = 0): 

d 2 k , a)/I . f 32 dk 2 lk 2 j . 

■Ar'(k)- / — n— = 0, Fla 



(2tt) 2 u v ' J 4vr 2nk 2 32 

d 2 ^ r2 Aa) /7 , f 32 dk 2 1 / k 2 , 

k 2 A[ a \k)+ ln--l =-— , (Fib) 



(2tt) 2 1 w i 4tt ttP V 32 / 4tt 2 ' 



(2tt) 2 p 7o 4tt 27rfc 2 V 32 y 2tt 

d 2 fc f 32 dk 2 1 fk 2 \ 1 , „ 

4 e) W-/ + l =--. (Fid) 



(2?r) 2 1 w 7 4vr vrA; 2 V 32 / 4vr 2 

At the same order of approximation, a number of intrinsically finite lattice integrals must 
be computed. When C2 = 0, some integrals can be evaluated analytically: 

(j2/>; k 2 A^\k) = - 1 -, (F2a) 



(2tt) 2 u v ' 4 
d 2 A; (fl) „, 3 



(2jr)a ^(*) = -4- ^ 

However some computations can only be, to the best of our knowledge, performed numeri- 
cally Hfl : 



G [ f ] = ~\j = 0.04616363, (F3a) 



^--J^^^k^O^,- (F3b) 
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FIGURES 

FIG. 1. Renormalization-group trajectories in the {/3 g ,/3 v } plane. 
FIG. 2. The large- iV mass ratio m^/mo, plotted as a function of x. 
FIG. 3. Feynman rules for the 1/N expansion in the continuum. 

FIG. 4. The 0(1) (subleading) finite scaling part of the free energy cp, computed from Eq. ( |5.4j) 
(solid line) and from Eq. (|5.6| ) (dot-dashed line); the 0(1/N) (subleading) finite part of the mass 
gap c m , computed from Eq. ( |5.24| ) (dashed line) and from Eq. ( |5.26| ) (dotted line). 

FIG. 5. 0(1/N) contributions to the two-point function. 

FIG. 6. The 0(1/ N) contribution to 5m 2 R , as a function of x. 

FIG. 7. Contributions to , the correlation function of the composite operator Pij(x). 

FIG. 8. Effective vertices in the continuum. All momenta are entering in the diagrams. 

FIG. 9. cz, the finite part of the renormalization constant Zp. 

FIG. 10. The static potential V(R), for several values of x. 

FIG. 11. 0(1/N 2 ) contributions to the topological susceptibility. 

FIG. 12. Feynman rules for the 1/N expansion on the lattice. 

FIG. 13. and a\ as functions of C2 («i is multiplied by 10 3 ). 

FIG. 14. Contour plots of A { Q a) (k) and A { 9 \k) + 2/vr. Contour lines are in logarithmic scale, 
separated by a factor \[2. 

FIG. 15. /io(0) (solid line) and Hq(\tt), (dashed line), normalized to thq. 
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FIG. 16. 0{\/N) contributions to the self-energy. 



FIG. 17. The renormalization-group flow trajectories of the lattice models, as defined by 



Eq. (11.28). 



FIG. 18. 8m\ as a function of /, for C2 = and k = 0.01, k = 0.1, k = 1/tt, and k = 1. Solid 
and dashed lines are the results of Eq. ( |12.1| ) for 9 = and 9 = 7r/4 respectively; dot-dashed lines 
are the results of Eq. ( 12. 5| ); dotted lines are the results of a power series expansion of Eq. (HJD 



to 3rd and 6th order. 

FIG. 19. bZ\ as a function of /, for C2 = and k = 0.01, k = 0.1, k = 1/tt, and k = 1. 
Solid lines are the results of Eq. ( 12.6j ); dotted lines are the results of a power series expansion of 



Eq. ( |12.6j ) to 3rd and 6th order. 

FIG. 20. 5m\ as a function of /, for the 0(2iV) model. Solid and dashed lines are the results 
of Eq. ( |12.4| ) for 9 = and 9 = 7r/4 respectively; dot-dashed lines are the results of Eq. ( |11.65| a); 
dotted lines are the results of a power series expansion of Eq. (|11.65| a) to 3rd and 6th order. 



FIG. 21. 8m\ as a function of /, for the CP N 1 model. The solid line is the result Eq. Qll.63 ); 



dotted lines are the results of a power series expansion of Eq. ( |11.63| ) to 3rd and 6th order. 

FIG. 22. The leading contribution to the difference Xt ~ Xt- 
FIG. 23. The tadpole graph contributing to the difference between Eq. ( |14.16| ) and Eq. ( |14.17| ). 

FIG. 24. The graphical form of the Schwinger-Dyson equations truncated to 0(1/N 2 ). Open 
circles indicate dressed propagators; full circles indicate dressed vertices 

FIG. 25. The rescaled energy levels r\ n as functions of t (solid lines). The dashed line corre- 
sponds to r) = t, where the excited levels disappear. 
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